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0.   Background 

This  report  is  one  of  a  series  of  studies  stemming  from 
current  research  on  electrohydrodynamic  (EHD)  power  genera- 
tion, sponsored  by  the  Department  of  Energy.   A  certain  EHD 
scheme  proposed  by  the  Marks  Polarized  Corporation  is  presently 
under  development.   The  authors  of  this  report  have  partici- 
pated as  consultants  on  various  theoretical  and  analytical 
aspects  of  this  effort.   One  of  our  main  objectives  has  been 
to  make  an  independent  appraisal  of  the  performance  possi- 
bilities and  limitations  inherent  in  the  general  type  of  EHD 
system  involved  in  this  development.   Another  objective  has 
been  to  create  a  sound  analytical  basis  for  optimizing  the  key 
design  parameters  of  such  a  device,  and  to  suggest  design 
improvements  where  appropriate. 

The  results  of  our  first  study  are  reported  in  detail  in 
Ref.  (1)  and  are  summarized  more  concisely  in  Ref.  (2).   This 
work  indicated  that  the  design  as  initially  proposed  would 
yield  unacceptably  poor  performance  and  suggested  a  certain 
design  revision.   Further  work  was  done  to  analyze  such  a 
revised  design  and  is  reported  in  the  appendix  of  Ref.  (1). 
An  improved,  simplified  and  shortened  version  of  this  work  is 
reported  in  Ref.  (3).   Unfortunately,  our  results  indicated 
that  the  revised  design,  wnile  showing  somewhat  improved  per- 
formance, still  would  not  attain  competitive  overall  efficiency 

The  Marks  Polarized  Corporation  nas  challenged  these 
pessimistic  conclusions.   They  argue  that  in  our  analysis 


the  compressibility  of  the  medium  was  not  taken  into  account 
and  that  if  it  had  been,  the  predicted  performance  would  have 
been  drastically  improved.   We  have  good  reasons  to  conclude 
that  Marks'  argument  about  compressibility  is  fallacious.   In 
the  first  place  our  analysis  did  consider  all  of  the  major 
effects  of  compressibility;  only  certain  minor  effects  in  one 
part  of  the  system  were  neglected  in  order  to  simplify  the 
calculations.   While  a  more  elaborate  and  accurate  treatment 
of  these  neglected  effects  is  certainly  possible,  we  maintain 
that  this  refinement  will  not  suffice  to  change  the  generally 
pessimistic  conclusions  established  in  our  original  analysis. 

Fortunately,  this  difference  of  opinion  between  ourselves 
and  the  Marks  Polarized  Corporation  pertains  to  a  question  of 
fact  about  which  it  is  entirely  possible  to  obtain  objective 
evidence.   To  this  end  we  present  in  this  report  a  revised 
and  improved  analysis  which  takes  all  compressibility  effects 
exhaustively  and  minutely  into  account.   Calculations  based  on 
this  revised  analysis  should  eventually  settle  the  above 
question  one  way  or  the  other.   A  typical  trial  calculation  is 
presented  later  in  this  report.   This  initial  calculation  has 
been  made  manually,  but  systematic  further  work  will  be  done  by 
computer.   The  full  results  and  conclusions  derived  from  these 
calculations  will  be  reported  later. 

It  has  also  been  conjectured  by  Marks  and  others  that  it 
might  be  possible  to  improve  greatly  the  performance  of  the 


EHD  generator  by  employing  a  two-fluid  cycle.   We  therefore 
present  in  this  report  a  detailed  analysis  and  procedure  by 
means  of  which  the  performance  of  any  specified  two-fluid  EHD 
system  can  be  calculated.   Our  analysis  takes  into  account 
fully  the  compressibility  of  both  fluids  at  all  locations  in 
the  system. 

We  propose  to  undertake  systemic  computer  calculations 
of  this  kind  which  should  enable  an  objective  assessment  to 
be  made  of  the  performance  potentialities  and  limitations  of 
several  proposed  two-fluid  EHD  generators.   Such  results  would 
be  of  great  value  in  helping  decide  whether  the  two-fluid  EHD 
generator  shows  sufficient  promise  to  warrant  further  develop- 
ment effort  and,  if  so,  what  the  optimum  design  point  values 
of  the  key  parameters  should  be. 

As  far  as  we  know,  the  present  report  is  the  only  pub- 
lished analysis  of  the  two-fluid  EHD  generator  which  includes 
the  compressibility  of  both  fluids  and  which  considers  the 
complete  thermodynamic  cycle.   The  principal  earlier  effort 
along  tnese  lines  is  given  in  Ref.  (4) ,  but  the  present  report 
goes  well  beyond  anything  undertaken  in  Ref .  (4) .   In  particu- 
lar, Ref.  (4)  makes  no  attempt  to  analyze  the  complete  thermo- 
dynamic cycle.   Moreover,  although  it  contains  certain  useful 
data  and  is  correct  in  some  respects,  Ref.  (4)  is  flawed  by  a 
number  of  serious  errors. 


While  the  present  report  is  otherwise  self-contained,  it 
does  assume  that  the  reader  has  a  general  familiarity  with  the 
type  of  EHD  power  generator  described  in  Refs.  (1) ,  (2)  and 
(3).   It  is  recommended  that  the  reader  unfamiliar  with  this 
general  scheme  first  review  briefly  the  introductory  material 
in  the  first  section  of  any  one  of  these  references  before 
proceeding  with  the  more  detailed  treatment  in  this  report  as 
this  earlier  background  information  is  not  covered  or  repeated 
in  the  present  report.   We  do,  however,  include  here  a  schema- 
tic diagram,  Fig.  0.1,  which  is  taken  from  Ref.  (3).   It 
illustrates  one  particular  version  of  an  EHD  power  generator. 
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1.    Introduction 

For  purposes  of  analysis,  the  overall  thermodynamic  pro- 
cess in  the  two-fluid  EHD  generator  may  be  represented  by  the 
schematic  flow  diagram  shown  in  Fig.  1.1. 

Primary  fluid  enters  the  ejector  at  station  1,  secondary 
fluid  enters  at  station  2  and  the  mixture  leaves  at  station  3. 

The  mixture  then  passes  through  a  diffuser,  process  3-*4, 
the  basic  purpose  of  which  is  to  reduce  the  kinetic  energy  of 
the  flow  and  thereby  decrease  the  associated  friction  losses 
downstream  of  the  diffuser.   On  the  other  hand  the  diffuser 
introduces  certain  losses  of  its  own  so  that  some  care  is 
necessary  to  optimize  the  design  in  this  respect. 

The  gas  mixture  then  passes  through  the  electrical  power 
section,  process  4+5,  from  which  the  gross  electrical  power  is 
extracted. 

Finally,  the  gas  mixture  enters  the  condenser/separator 
at  station  5.   The  primary  fluid  is  separated  from  the  gas 
mixture  by  condensation  and  leaves  at  station  6.   The  fluid  at 
this  point  is  assumed  to  be  compressed  liquid  at  a  known  static 
pressure  Pg   and  a  known  static  temperature   Tfi  .   Kinetic 
energy  at  station  6  is  regarded  as  negligible.   The  secondary 
fluid  leaves  the  condenser/separator  and  enters  the  ejector  at 
station  2.   We  assume,  moreover,  that 


and 


P6  =  P2  =  Px  (1.1) 


T6  =  T2  (1.2) 
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The  primary  fluid  which  leaves  the  condenser/separator 
at  station  6  is  circulated  back  through  the  pump  and  through 
the  boiler/superheater,  process  7+1. 

It  should  also  be  explained  that  the  primary  fluid  at 
station  1  is  a  condensible  vapor  of  high  molecular  weight 
which  contains  fine  liquid  droplets.   Unless  specifically 
stated  otherwise  the  quality   z   at  this  location  is  taken 
as  0.95.   The  important  thermodynamic  static  properties  at 
station  1  can  be  expressed  in  the  following  form. 

h,  =  hf  +  z  hf    Enthalpy  (1.3) 

s,  =  s.p  +  z  sf    Entropy  (1.4) 

Here  subscript   f   denotes  the  saturated  liquid  and  subscript 
fg  denotes  the  change  between  saturated  liquid  and  saturated 
vapor.   All  properties  are  evaluated  at  a  specified  pressure 
P,  ;  the  corresponding  saturation  temperature   T,   is  then 
also  known. 

The  corresponding  stagnation  properties   T  ,   and  h  , 
are  fixed  when  the  stagnation  pressure   P  -,   is  specified. 
Of  course  the  entropy  has  the  same  value   s,   at  the  stagna- 
tion condition  as  at  the  corresponding  static  state.   Thus 

specification  of  P  ,   and  s,   suffices  to  fix  all  other 
e  si        1 

stagnation  properties. 

While  it  is  convenient  to  start  the  analyses  by  stipula- 
ting the  quality  at  station  1,  it  is  also  advisable  later  in 
the  calculation  to  determine  the  corresponding  quality  at 


station  4,  the  entrance  to  the  electrical  working  section. 
It  is  essential  that  the  fluid  entering  the  working  section 
contain  finely  dispersed  liquid  droplets  of  the  proper  size. 
These  droplets  carry  the  electrical  charges  that  are  the 
central  feature  of  EHD  power  generation.  We  assume  tentatively 
that  acceptable  droplet  characteristics  are  obtained  by 
maintaining  the  quality  at  the  inlet  to  the  working  section 
in  the  range  0.92  to  0.98.   These  limits  are  only  estimates. 
The  thermodynamic  analysis  would  be  greatly  simplified 
if  both  fluids  were  perfect  gases  with  constant  specific 
heats.   In  fact  only  the  secondary  fluid  satisfies  this  re- 
quirement.  Nevertheless,  we  can  define  an  "equivalent  perfect 
gas"  which  adequately  approximates  the  essential  thermodynamic 
properties  of  the  primary  vapor/liquid  mixture.   When  carried 

out  judiciously,  this  procedure  simplifies  the  analysis  at  the 

5 
cost  of  only  a  small  loss  in  numerical  accuracy  . 

The  ratio  of  specific  heats   y   of  the  "equivalent  perfect 

gas"  may  be  deduced  from  the  important  relation 

Yi/CYi  ~  1) 
(Pal/Pl)  =  (Tsl/Tl}  (1'5) 

Solving   this   for      y,      gives 

,-1 


Yl  = 


_    *n(Tsl/Tl> 


*n(Psl/Pl) 


(1.6) 


Another  fundamental  property  of  the  "equivalent  perfect 
gas"  is  its  specific  heat   C  ,   which  may  be  evaluated  from 
the  expression 


<hsl  -  hi> 

Cpl    (Tgl  -  Tx)  (1*7) 

where  all  quantities  on  the  right  are  now  known. 

It  then  follows  from  standard  perfect  gas  relations  that 
the  gas  constant   R,   and  the  molecular  weight  W,   of  the 
"equivalent  perfect  gas"  are,  respectively, 

"l-ftrrv'ri  (1-3) 

and 

w  =  I  (1-9) 

1    Ri 

where   R  =  universal  gas  constant 
=  8315  joules/kg  °K 

Eqs.  (1.6)  through  (1.9)  complete  the  definition  of  the 
"equivalent  perfect  gas"  which  adequately  simulates  the  essen- 
tial thermodynamic  properties  of  the  primary  vapor/liquid 
mixture. 

It  should  also  be  explained  that  the  secondary  fluid  at 
station  2  is  predominantly  a  noncondensing  gas  of  low  molecu- 
lar weight  primary  vapor. 

It  can  be  shown  that  the  mass  ratio  oj  of  condensible 
primary  vapor  to  dry  secondary  gas  at  station  2  is  given  by 
the  expression 

Wv    W 

"=  ^ip2- wj  (1-1) 
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where 

W  =  true  molecular  weight  of  condensible  vapor  (not  to 
be  confused  with  the  "equivalent  molecular  weight 
W-,"  considered  earlier) 

W„  =  molecular  weight  of  noncondensing  gas 
P  (Tfi)  =  vapor  pressure  of  condensible  primary  fluid  at  known 
temperature   TV 

P~  =  known  static  pressure  at  station  2 

The  present  analysis  is  restricted  to  conditions  under 
which   co   is  very  small  compared  with  unity  and  may  be  neglec- 
ted.  Thus  the  fluid  at  station  2  may  be  treated  as  dry  gas. 
This  greatly  simplifies  the  analysis.   Eq.  (1.10)  is  useful  for 
verifying  that  this  assumption  is  indeed  satisfied  in  any  par- 
ticular instance. 

It  should  be  added  that  if  it  later  becomes  of  interest 
to  make  calculations  for  circumstances  under  which  parameter 
oo   is  not  negligible,  the  present  analysis  can  in  fact  be 
generalized  to  include  this  effect.   For  the  present,  however, 
we  prefer  to  deal  only  with  the  simpler  situation  in  which  co 
may  be  neglected. 
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2.    Ejector 

The  ejector  receives  a  primary  stream  of  high  molecular 

weight  gas  at  station  1,  a  secondary  stream  of  low  molecular 

weight  gas  at  station  2  and  discharges  the  resulting  mixture 

at  station  3.   Static  conditions  at  these  three  stations  are 

designated  by  subscripts   1  ,   2  ,   3  .   Stagnation  conditions 

at  the  corresponding  stations  are  designated  by  subscripts 

si  ,   s2   and  s3  . 

The  following  quantities  are  arbitrarily  specified  or 

known  at  stations  1  and/or  2,  namely, 

Molecular  weights:   W,   and  W2 

Ratios  of  specific  heats:   y-i   and  y2 

Stagnation  pressure  of  primary  jet:   P  , 

Stagnation  temperature  of  primary  jet:   T  , 

Static  pressures:   P,  =  P2  (=  Pg) 

Static  temperature  of  primary  jet:   T,  (=  saturation  tempera- 
ture of  primary  fluid  at  pressure  P,) 

Static  temperature  of  secondary  jet:   T~  (=  Tfi) 

Mass  flow  ratio:   m,/m~  =  m,/(m,  +  nu)  =  x 

Velocity  ratio:   v-]/vo  =  ^ 

We  treat  the  fluids  at  stations  1,  2,  3,  4  and  5  as  per- 
fect gases  with  constant  specific  heats. 

The  velocity  ratio  of  the  ejector  may  be  developed  as 

follows: 


li  -  v  -  !i  ^  JTi  RiTsi      Mi       _i  ,,  ,, 

V2     y     a2   M2  -W  T2  R2  T2  A|Tsl/Tl  •  M2 
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We  can  now  use  this  result  and  other  standard  perfect 
gas  relations  to  complete  the  calculation  of  key  thermody- 
namic properties  at  stations  1  and  2.  It  is  convenient  to 
arrange  these  calculations  sequentially  as  follows : 


Mi  =  VTTT^-rr  (*r ' ')  (2-2) 


M2  =  ( V 


$ 


II  Rl  Tsl    Ml 
Y2  R2  T2  U  Tsl/Tl 


(2.3) 


^«s2  "^2  "       2 

ff  '   1+   "iT-  M2  (2-4' 

Ts2  "  T2  (Ts2/T2)  (2-5) 


(2.6) 


PS2  =  P2  (gf)  (2-7» 

Before  we  can  determine  the  actual  conditions  at  station 
3,  the  outlet  of  the  ejector,  it  is  first  necessary  to  analyze 
two  other  hypothetical  cases  as  indicated  schematically  in 
Fig.  2.1.   In  all  three  of  these  cases  the  flow  is  taken  as 
steady  and  adiabatic. 

The  first  case,  that  shown  in  Fig.  2.1(a)  ,  involves  a 
hypothetical  device  which  receives  two  gas  streams  at  stations 
1  and  2  and  discharges  the  resulting  mixture  at  station   x  .   In- 
let conditions  at  1  and  2  are  identical  to  those  of  the  actual 
ejector. 


13 


w    (l-x) 


© 


Hypothetical  Minimum 
Loss  Device 


(a) 


© 


© 


^ 


(?)  ** 

V    (l-x) 


© 


Actual   Ejector 


(c) 


© 


Fig.    2.1    COMPARISON  OF    ACTUAL   EJECTOR   WITH   IDEAL 
EJECTOR    AND  WITH  HYPOTHETICAL  MINIMUM 
LOSS  DEVICE. 
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Recall  that  the  adiabatic  mixing  of  two  different  gases  is 
an  inherently  irreversible  process  that  always  involves  a  cor- 
responding entropy  increase.   We  stipulate  that  the  first  case 
is  such  that  the  only  irreversibility  which  occurs  is  that  asso- 
ciated with  this  mixing.   We  term  this  mass  mixing  to  distin- 
guish it  from  another  type  of  mixing  considered  below  which  we 
term  momentum  mixing. 

The  second  case,  that  shown  in  Fig.  2.1(b) ,  involves  an 
ideal  ejector  which  receives  two  gas  streams  at  stations  1  and 
2  and  discharges  the  resulting  mixture  at  station   y.   Again, 
the  inlet  conditions  at  1  and  2  are  identical  to  those  of  the 
actual  ejector.   This  case  satisfies,  among  other  relations,  the 
idealized  one-dimensional  momentum  equation  for  frictionless , 
constant  area  flow.   It  can  be  shown  that  in  this  case  there  is 
an  overall  entropy  increase  which  includes  not  only  the  previ- 
ously mentioned  effect  of  mass  mixing  but  also  a  further  increase 
associated  with  momentum  mixing. 

The  third  case,  that  shown  in  Fig.  2.1(c)  ,  represents  the 
actual  ejector  itself  which  receives  input  streams  at  stations 
1  and  2  which  discharges  the  resulting  mixture  at  station  3. 

All  three  of  the  above  cases  refer  to  constant  area  pro- 
cesses in  the  sense  that 

Ax  =  Ay  =  A3  =  (A1    +    A2)  (2.8) 

Once  the  hypothetical  operating  conditions  at  stations  x 
and  y  have  been  found,  it  then  becomes  possible  to  express  the 


15 


corresponding  actual  conditions  at  station  3.   This  is  accom- 
plished through  the  use  of  an  ejector  effectiveness   nF  . 
This  parameter  is  defined  more  explicitly  in  the  later  analysis. 
Values  of  x)„     must  be  estimated  by  reference  to  test  data  on 
ejectors  (i.e.,  Ref.  6,  7,  8). 
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3.   Adiabatic  Mixing  of  Two  Streams 

The  process  in  the  ejector  is  treated  as  the  adiabatic 
mixing  of  two  streams.   The  following  relations  can  be  shown 
to  apply. 

R  =  8315   Joule/kg°K  (3.1) 

R^f-  (3.2) 

R2=|-  (3.3) 

R   =  R   =  R.  =  x  R..  +  (1  -  x)  R.  (3.4) 

X     y     3       1  2 


(3.5) 


cpl  = 

Yl 

Rl 

<*1 

-    1) 

CP2   = 

Y2 

R2 

(Y2 

-    1) 

(3.6) 


C=C   =C,  =  xC,  +  (l-x)Cn  (3.7) 

px     py     p3       pi  p2 

c  3 

^x  -  ^y  =  ^3  =  (Cp3P-  R3)  (3'8) 

Tsx  =  Tsy  "  Ts3  -  C^  [X  CP1  Tsl  +  U  r  X)  CP2  Ts2]    ^'^ 
Ax  =  A   =  A3  =  (Ax  +  A2)  (3.10) 

Notice  that  the  quantities   R  ,  C   ,  y    ,    T    and  A   all 

*  p  s 

have    the   same    values    at   stations    x,    y   and   3. 

This    fact  may  be   used   to   simplify   many   of   the   subsequent 
expressions.      It  is    always   permissible   to   substitute   subscript 
3    for   subscripts    x  or  y   on    any   of    these   quantities. 
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4.      Mass    Flow 

It   is    convenient   temporarily    to   omit   station   subscripts 
1,    2,    3,    x,    y   and   thereby   develop    certain   needed  mass    flow 
relations    in   generalized   form.      Thus    the   mass    flow   across    an 
arbitrary   station   may  be  written    and   developed  as    follows. 

m   =    pAV  =  ( |^j  A    (/7RT   m)    =    PAM  ^ 


P    .A 
s 


T    /T 
s 


Vrttt  <Vp> 


M 


(4.1) 


This   may  be    rewritten   as 


P      A 
m  =      — - f  (M) 


JM^y 


(4.2) 


where    the    auxiliary    function      f(M)       is    defined   as    follows. 


It  /t 
f(M)  =  V7PT  M 


=        M 


i  +  (1^-i)  m^ 


(Y    +    1) 
2(Y    -    1) 


(4.3) 


Eq.  (4.2)  can  be  applied  specifically  to  stations  1,  2, 
x,  y,  3  as  follows 


P   A 
m,  =    S1  L     ■   f(Mn) 

'  1  sl/yl 


(4.4) 


P  ,  A 
m,  =  -— ^i— ±-  ■  f(M  ) 

2   s2'  '  2 


(4.5) 
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P         A- 
K  =  SX  f (M)  (4.6) 

^R3    Ts3^3 


P         A 
m     =      ^ —   f(M   )  (4.7) 

"  V*3   Ts3^3  y 


P    -    A 
A3   =  S         \         f(M3)  (4.8) 


VR3   Ts3^3 


also 


m 

m~l  -  x  (4.9) 

m3 

±     =  ±     =  *      =  m     +  ft.  (4.10) 

x  y  3  1  2 

From  Eqs.     (4.4),     (4.5),     (4.9)    and    (4.10)    we   readily    find 
that 


'tk\  -  x  /PS2\  /Rl   Tsl   Y2    f(M2> 

A2 )  '     TT^1^   \Psl//R2    Ts2    ~1   W 


Then    from  Eq.     (2.10)    we   infer  also    that 


(4.11) 


\\  (A^) 


A3/        "(1   +    A1/A2) 


(4.12) 


Since   all   quantities   on    the  right  side   of  Eq .     (4.11)    are 
now  known,    Eqs.     (4.11)    and    (4.12)    fix  the    area   ratios      at/A2 
and     A  /A_    . 

From  Eqs.     (4.5),     (4.6),     (4.9),     (4.10)    and    (2.8)    we    deduce 
further   that 
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(1  -  A  /A  )  /P   \   R   T  -  y0 

«V  "  ^r^  (fg)  J*^ f  "*       (4-13) 

All  quantities  on  the  right  side  of  Eq.  (4.13)  except 

P    are  now  known.   The  method  of  determining  P    is  ex- 
sx  3    sx 

plained  in  a  later  section;  refer  to  Eq.  (5.15).   Once   P 

S  X 

has  been  found,  Eq .  (4.13)  fixes   f(M  )   whereupon  Eq.  (4.3) 
fixes   M   itself.   The  solution  of  Eq.  (4.3)  for  M  when 
f(M)   is  known  involves  an  iterative  procedure  which  is  ex- 
plained later  in  this  section. 

From  Eqs.  (4.6),  (4.7),  (4.8)  and  (4.10)  we  also  find  that 

P    f(M  )  =  P    f(M  )  =  P  ,  f(M_)  (4.14) 

sx    x     sy    y     s3    3 

Assuming  P    and  M   known,  this  relation  along  with  Eq.  (4.3) 

SX  X 

fixes   P    when   M   is  specified;  it  also  fixes   f(M_)   and 
sy         y  3 

M_   when   P  ->   is  specified.   The  methods  of  calculating  M 
3         s3      r  y 

and  P  -.   are  explained  later;  refer  to  Eqs.  (6.8)  and  (7.14)  . 

Fig.  4.1  is  a  rough  sketch  of  the  function  f(M)  as  defined 
by  Eq.  (4.3)  .  This  function  vanishes  at  M  =  0  and  at  M  =  «  . 
The  maximum  value  occurs  at   M  =  1   and  equals 

(Y  +  1) 

,.|2 
f 


max 


Y  +  1 


2(Y  ~  1) 

(4.15) 


It  is  clear  from  the  figure  that  in  the  range   0<f<  f  „.,  , 

3  max 

Eq.  (4.3)  has  two  positive  real  roots,  one  subsonic  and  the 
other  supersonic.   Usually  the  supersonic  root  must  be  discarded 
because  it  is  not  consistent  with  a  further  constraint  imposed 
by  the  second  law  of  thermodynamics. 
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t 

f(M) 


M  > 


Fig.  4.1  APPROXIMATE  SKETCH  OF  FUNCTION  f(M) 
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In  order   to   solve   Eq.     (4.3)    for     M     when   values   of      f 

and     y      are   specified,    we   utilize    an   iterative   procedure  based 

on  Newton's   method.      Let     M        denote    the   nth    trial   value   of 

n 

M    .      We   then   define    a    function      F(M   )       and  its    derivative 

n 

i 

F  (M  )   as  follows, 
n 


F(M  ) 
n 


=  f 


1  + 


Y  -  1 


M 


n 


(Y  +  l) 
2(y  -  l) 


-  M 


n 


(4.16) 


F*  (M  ) 

n 


f(i^ 


M 


n 


1  + 


Y  -  1 


M 


n 


(3  - 


2TT= 


Y) 

T7 


-  l 


(4.17) 


The   (n  +  l)3t  trial  value  of  M  may  now  be  taken  as 


F(M  ) 

M,     ,v  =  M„  -  ■  !  n 

(n  +  1)     n 


(4.18) 


F  (M  ) 
n 


The    cycle   of   calculations    defined  by  Eqs .     (4.16)  ,    (4.17) 
and    (4.18)    is    repeated   until    the    result   for     M        converges    to 
a   stable   value    at   the    desired   level   of   accuracy. 
Thus 


M  =    Him  M 


n 


(4.19) 


It  is    readily   apparent      that   the   value  of     M     so   obtained 
satisfies   Eq.     (4.3)     for   the  prescribed  values   of      f      and     y    . 
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5 .   Entropy 

Consider  the  case  shown  in  Fig.  2.1(a).  The  specific 
entropy  s  of  the  discharged  gas  may  be  expressed  in  the 
form 


sx  =  x 


*pl 


X 


■  r.  I  |5 


xl 


n 


+  (1  -  x) 


C  0  in 
p2 


2      n 


x2 


(5.1) 


where   P  ,   and  P  0      are  the  partial  pressures  of  the  two  com- 


xl 


x2 


ponents   of   the   gas   mixture.      Symbols      P        and      T        denote    the 
pressure    and  temperature   of   the   ambient  atmosphere.      The   en- 
tropy of  each   component  is    assigned   the  value    zero    at   this 

reference   state      P      ,    T 

o  o 

Incidentally,  the  primary  fluid,  being  condensible,  might 
well  exist  only  in  the  liquid  state  at  the  reference  condition 
P      ,    T      .      Nevertheless,    so   long  as    it   is    in    the   gaseous    state 

O  O  '  z>  z> 

at  station     x    ,    Eq.     (5.1)    may   still  be    used. 

The  partial  pressures    of    the  two    components    are   propor- 
tional   to    the   respective   mol    fractions.      Hence 


xl 


X    R1 

X    R      +     (1    -    x)R2 


X    R.    /P 
_1    /      X 


(5.2) 


x2 


(1    -   x)R2 
x   Rx   +    (1   -   x)R2 


(1   -   x)R, 

rT 


X 


(5.3) 
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Upon   substituting  Eqs.     (5.2)    and    (5.3)    into    (5.1)    and 
simplifying  we   may  obtain    the   result  in   the    form 

sx  ■    sm  +  Cp3   £n  (¥)  -   R3   £n  (**)  (5-4» 

where 
s^  =   R      In   R3   -   x   Rx    £n    (x   ^)    -    (    1   -   x)R2    £n    [(1   -  x)R2] 

(5.5) 

It  can  be  shown  that  this  quantity   s    represents  the 

m    c 

entropy  increase  associated  with  the  mass  mixing  of  the  two 
different  gases. 

Moreover,  since  the  entropy  of  the  stagnation  state   s 
is  by  definition  identical  to  that  of  the  corresponding  static 
state   x  ,  we  may  replace   T   and  P    in  Eq.  (5.4)  by   T 
and  P    ,  respectively.   Thus  we  obtain 

sx  =  sm  +  Cp3  tn  (ff)  "  R3  £n  (If)         (5-6) 

In  previous  studies  of  this  series,  we  have  considered 
only  the  special  case  in  which  the  same  fluid  is  employed  for 
both  the  primary  and  secondary  streams.   In  that  case  we  have 

R±   =   R2  =  R   =  R  (5.7) 

C,=C0=C,  =  C  (5.8) 

pi     p2     p3     p 

y1   =  Y2  =  Y3  =  Y  (5.9) 

Moreover,  in  this  special  case  there  is  no  mass  mixing,  so 
that  we  must  set 
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s   =  0  (5.10) 

m 


Consequently,  Eq.  (5.6)  now  simplifies  to 


sx  =  o+  cP  *n(H  -  R£n(iH      <5-ii) 


Next  reverting  to  the  case  of  two  different  gases,  let 
At  denote  a  ,time  interval  during  which  unit  mass  crosses 
station  x  .   During  this  same  interval   x   units  of  mass  cross 
station  1  and   (1  -  x)   units  of  mass  cross  station  2.   The 
entropy  which  leaves  across  station   x  must  equal  the  entropy 
which  enters  across  stations  1  and  2  plus  the  entropy  increase 
caused  by  the  mass  mixing.   Thus 

s   =s   +  x  s,  +  (1  -  x)s~  (5.12) 

x    m      1  2 

where 

si =  V  ln  C^r)  '  Ri  *n  l^1)  (5-l3) 


s2  "  Cp2  £n  (ff)  "  R2  ln   (Jf)  (5-14) 

Upon  eliminating  s    between  Eqs .  (5.6)  and  (5.12),  we 

A 

find   that     s        also    cancels    from  the   result.      It  is    then   a 
m 

simple   matter  to   solve    for   the    exit  stagnation  pressure   in    the 
form 

*n  (H  =  k  icP3  *n  Gn  -  *■  si  - (1  -  x)s2}    <5-i5) 
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Thus  the  sequence  defined  by  Eqs .  (5.13),  (5.14)  and 
(5.15)  now  fixes   P    .   Next  referring  back  to  Eq.  (4.13), 
we  can  calculate  the  value  of   f  (M  )  ;  finally,  from  Eqs. 

A 

(4.16),  (4.17),  (4.18)  we  can  calculate   M   itself.   Once 

ft 

P    ,  T   =  T  0   and  M   are  known,  it  is  a  simple  matter  to 
sx    sx    s3        x  r 

calculate  the  corresponding  conditions   P   ,  T    from  the 

X       X 

usual   perfect   gas    relations. 
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6.      Momentum 

Consider   the   idealized  one   dimensional    constant  area 
ejector      shown   in   Fig.    2.1(b).      The   momentum  equation    for   this 
system  may  be  written 

(pl  +  pl  V12)A1  +   (p2  +  p2  V22)A2  =   (Py  +  py  Vy2)Ay      (6>1) 

The   mass    conservation    relation    for    this    system  can  be 
expressed  in   the   form 

Pl   h     \   .    »2    *2    V2    .    p      A      v  (6.2) 

x  (i  -  x)         y    y    y 

Let  us  divide  the  three  terms  of  Eq.  (6.1)  by  the  corres- 
ponding three  terms  of  Eq.  (6.2).   Notice  that  the  areas  can- 
cel out  of  the  result.   In  this  way  we  obtain 

(Pl  +  Pl  v,2)  (p2  +  p2  v22)  _  (Py  +  py  vy2) 

Pl  Vl        (      }      P2  V2      "      ^  Vy 

(6.3) 
Temporarily   dropping   the   station   subscript,    we   next  deve- 
lop  the   typical    term  of  Eq.     (6.3)    in   generalized   format  as 
follows. 

(P    ±    pV2)     u   [P    *      h       (^RT)M2]    =     FT   (1    +   YM2) 

pV  L,    /Trt  m  v  Y      m  /t-7t 


RT  '  s' 


JRT 


RTS  (1    ±    YM2)  (6'4) 


Observe   that    the  pressure      P      cancels    from  the   result. 
It  is    convenient   to   define   the   auxiliary   function 
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g(M)  - 5-  A    +  (3L-I_Jl)M2  (6.5) 

(1  +  yM  ) 

With  this  notation  Eq.  (6.3)  may  be  rewritten  in  the 
form 


x   /Rl  Tsl    (1  -  x)   7*2  Ts2     1     7*3  Ts3 


glM]T/      yx    '  +    gTMp     /     y2  glM^T  /     y3  (6,6) 

Rearranging   gives 


7*3   Ts3  i        x  Al    T_sl         (1   -   x) 


g(V   "  /      y3       7TV    /      y1  "gTM^T 

Since    all   quantities   on    the    right  are   known,    Eq.     (6.7) 
fixes      g(M  )     .      Then     M        itself   can  be    found  by  inverting 

y  y 

Eq.  (6.5).   Fortunately,  an  explicit  solution  is  possible  in 

2 
this  case  as  Eq.  (6.5)  can  be  reduced  to  a  quadratic  in  M 

The  result  is 


,2  _  (1  -  2Yg2)  ±  /l  -  2(y  +  Dg: 
1  -  Yd  ~  2Yg2) 


M*  =  vx    M^  -yx    "vjj  Z   ±12  (6.8) 


The  general  character  of  the  function   g(M)   is  sketched 

in  Fig.  6.1.   Notice  that  g(M)   vanishes  at  M  =  0  ,  reaches 

its  peak  value   g     at  M  =  1  ,  and  decreases  toward  the 
v  ^max 

limit   g   as   M-*»  .   It  can  be  shown  that 


g    = (6.9) 

maX    ^2(Y  +  1) 


-\/^^  <«•"> 
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t 

g(M) 


M ► 


Fig.  6.1  APPROXIMATE  SKETCH  OF  FUNCTION  g(M) 
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Fig.  6.1  discloses  the  existence  of  a  supersonic  root  over 
the  range   g  <g<g  „  , .   This  root  should  not  be  accepted  until 
a  check  calculation  is  made  to  determine  whether  the  result  is 
consistent  with  the  second  law  of  thermodynamics.   The  subsonic 
root  is  found  by  retaining  only  the  negative  sign  before  the 
radical  in  Eq.  (6.8). 

Once   M   has  been  found  from  Eq.  (6.8)  ,  the  corresponding 

value  of  P    may  be  found  from  Eqs .  (4.3)  and  (4.14).   This 
sy    J 

suffices  to  fix  all  properties  at  states   y   and  sy  . 
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7.   Availability 

The  entropies   s1   and  s~   of  the  two  input  streams 
were  defined  earlier  in  Eqs.  (5.13)  and  (5.14).   The  corres- 
ponding steady  flow  availability  functions  with  respect  to  an 
ambient  atmosphere  at  pressure   P    and  temperature   T   may 
be  written 

*I  =  CP1  (Tsl  -  V  "  To  sl  (7-1) 

*2  "  Cp2  (Ts2  '  V  "  To  Sl  (7-2> 

Consequently   the    total   available   energy   entering   the   sys- 
tem becomes 

^m  =   x    ^1    +    (1    -   x)^2  (7.3) 

The    respective   entropies   of   the   streams    leaving  at  stations 
x    ,    y      and   3      are 

sx  ■   Sm  +  Cp3    ^  (ff)   "   R3    ln  (5f )  (7-4) 

sy  "   sm  +   Cp3    ln  (ff)  "    R3    to  (if )  ('7'5) 

s3  "   sm  +  Cp3    ln  (ff)  "    R3   *n  (if)  (7"5) 

where   s   has  been  previously  defined  in  Eq.  (5.5). 

The  corresponding  availabilities  of  the  streams  leaving 
at  stations  x,  y  and  3  may  therefore  be  written 
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tJ>      -CQ(T,-T)-T      s  (7.7) 

yx  p3        s3  o  ox  v    .    / 

f      =C,(T,-T)-T      s  (7.8) 

yy  p3        s3  o  o      y 

*,  =   C    -    (T    -   -   T   )    -   T      s.  (7.9) 

r3  p3        s3  o  o      3 

Now   consider   the   losses   of   availability   listed  below. 
By   substituting   the    above   expressions    for   the   quantities   on 
the    left  and  simplifying  we    readily   obtain   the   expressions 
shown   on    the    right.      Thus 

<*m  "   *x>    "  To   sm  (7-10> 

(*x"   V    "   R3   To    ln(^)  <7al) 

<*x"   *3>    "   R3  To    4n(?Jf)  <7,12) 

These    results    are   very   significant.      Eq.     (7.10)    defines 
the   loss   of   available   energy   caused  by   the  mass   mixing  of   the 
two   different  gases.      Eq.     (7.11)    defines    the    further   loss   of 
available   energy   caused  by    the   momentum  mixing   in   an   idealized 
ejector.      Eq.     (7.12)    defines    the    augmented   loss   of   available 
energy   caused  by    the   momentum  mixing  in    the   real   ejector. 

It  is    useful   to   postulate    an   empirical    relation  between 
these    last   two   losses  which   is    shown  below   on    the   left.      The 
expression   on    the    right   then    follows    from  Eqs.     (7.11)    and 
(7.12).      Thus 

(*x  -   *3)    -  ^E   -    ^    (Ps3/Psx)  <7-13' 
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We  term  nE   the  ejector  effectiveness.   Its  value, 
which  is  always  less  than  unity,  must  be  estimated  from  appro- 
priate test  data  on  ejectors. 

If  we  treat  nE   as  known,  Eq.  (7.13)  fixes   P  3  .   The 

solution  is  simply 

1_ 

nE 

P   =  P   [   sy_  ]  (7.14) 

s3    sx  \  P   / 
\  sx  / 

We   can  next   find      f(M_)       from  Eq.     (4.14)    and     M,      from 
Eqs.     (4.16),    (4.17)    and    (4.18).      This   suffices    to    fix  all  pro- 
perties   at  states    3    and  s3    thereby    completing,,  the   analysis   of 
the   ejector. 
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8.      Diffuser 

It  is   convenient   to    analyze    the   diffuser   in    terms    of   the 

concept  of   availability   in   steady    flow.      However,   we  wish    to 

deal  with    the    characteristics   of   the   diffuser   itself   and   these 

cannot  depend  on    the   arbitrary   values      P      ,    T       which  happen 

to    characterize    the   condition  of   the   ambient  atmosphere.      We 

can    achieve   our  objective   by   defining   the   availability  with 

respect   to    a   reference    state      P-.    ,    T-.     which   characterizes 

the    diffuser   rather    than  with    respect   to    the   ambient  state 

P      ,    T      . 
o    '      o 

Under   these    circumstances   we   may  write    the   appropriate 
availabilities    at  diffuser   inlet   and  outlet   as    follows 


*s3   "   Cp3    <Ts3  -    T3'    "    T3  lCp3    ln  (ff)-    R3    tn  (ff  )  }       (8-1) 


*s4    =  Cp3    (Ts3   -    T3>    "    T3      Cp3   £n  ft1)-    R3    ia(?r)l       <8-2) 


These   equations   make    use   of   the    fact   that     T    „    =   T    _ 
^  s4  s3 

Also,    states      s3      and      3      are    at   the    same   entropy   so   that 


<SS3   -    S3>    "    °    "   Cp    *"  (ff)   "    R3      ^  (Jr)  <8-3) 


From  Eqs.     (8.1)    and    (8.3)    we    conclude    that 

IP    ,    =    C    .     (T    .    -    T.)     =  (-4- J  (8.4) 

rs3  p3         s3  3  \     2     / 
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This    result   shows    that   the   available   energy   at   the 
diffuser  inlet   is    simply    the   inlet  kinetic  energy   itself. 

The   loss   of   availability   through    the   diffuser  may  be 
found  by   subtracting  Eq.     (8.2)     from    (8.1).      The   result  is 

(*s4   "    *s3>    -   A*  "   +   R3   T3    in(^j  (8-5) 

We   now   define   the  diffuser  effectiveness   as 

;H"nD  (8-6) 

S3/ 
so   that 

/Ps3 

-  ,         _    ,      \  R3   T3   in  f  —^ 

s3   "    ^s4\         Aip      _    n  »  \    s4/ 

T"; — J  -  j—  ~   (1  "  riD)    "  C  I    (T  I  -  TO  (8-7) 

rs3        /  rs3  p3         s3  3 


This    reduces    readily    to 

P 


s3 

In 


(i  -  n.J  = 


Ps4/  (8.8) 


D  /        »    2\ 


Solving  for  the  pressure  ratio  gives 

2 


Y3  M, 


Ps4       "  (1  "  nD)    "2- 
e 


Ps3  (8.9) 

This    is    the    result  required.      It   fixes      P    .      when      P    3    , 
M_      and     rin      are   specified. 
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Assuming   the   area   ratio     Ao/A4      to  be    arbitrarily   spec- 
ified,  we    can    find     M        from   the   continuity   relation.      Refer- 
ring  to   Eq.     (4.8)    we   may  write 

m-.   =  s3      3-    ■  f(M  )    =  s4      4     -  f(M   )    =   m  (8.10) 

^3   Ts3^3  ^3   Ts3'Y3 

Consequently 

f  ("4>  =  fe)  ft) f  <M3>  "-11' 

Eq.     (8.11)     fixes      f(M4)     •      Then      M.      follows    in    the    usual 
way    from  Eqs.     (4.16),    (4.17)    and    (4.18).      The    result  suffices 
to    fix  all   properties    at  states      4      and     s4    .      Thus 


W-GW-^P^W 


T4   =   Ts3/(TS3/T4>  (8-13) 


;« ■  &  ^ 


(8.14) 


P4    "   Ps4/(PS4/V  (8-15) 

The    configurations    studied   in   earlier  work   did  not   incor- 
porate   a   diffuser.      For   the  purpose   of   comparing   the    results 
of   the   present   analysis  with   corresponding  earlier  work,    it  is 
desirable   to  be    able    to  eliminate    the   effects   of   the    diffuser 
in   certain    cases.      This    can  be    accomplished   in    the  present 
analysis    simply  by   setting 
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nD  =  l  (8.16) 
and 

It  is   easy    to   confirm  that  under    these    circumstances 

Ts4  =  Ts3  (8.18) 

Ps4  =  Ps3  (8.19) 

M4  =  M3  (8.20) 

A4  =  A3  etc.                                           (8.21) 


and  all   diffuser  effects    disappear, 
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9.    Electrical  Power  Section 

The  power  conversion  process  in  an  EHD  duct  is  usually 
treated  on  the  basis  either  of  constant  area  or  of  constant 
static  state.  Inasmuch  as  electrical  power  output  per  unit 
mass  is  small,  the  numerical  differences  between  the  results 
computed  by  these  two  methods  is  negligible.  For  definite- 
ness  in  this  analysis,  however,  we  assume  constant  area. 

The  electrical  power  that  can  be  obtained  from  an  EHD 
duct  of  constant  area,  negligible  change  of  density  and  optimum 
length  can  be  estimated  from  the  one  dimensional  version  of 
Poisson's  equation  which  governs  the  electrical  field.   The 
solution  is  well  known  and  will  not  be  derived  here;  a  detailed 
derivation  may  be  found  in  Ref .  (3) .   The  essential  result  may 
be  written  in  the  form 

Pe=I£EbA4V4  {9-1] 


where 


P   =  gross  electric  power  output,  watts 

e   =  permittivity  of  medium 

-12 
=  8.854  x  10      farad/m  (for  any  gas) 

E,  =  dielectric  strength  of  medium,  volts/m 

A,  =  area  of  duct  (constant) 


V 


.   =  velocity 
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According  to  the  test  data  analyzed  in  Ref.  4,  the  dielec- 
tric strength  is  well  approximated  by  the  expression 

Eb  ■  CB3R3P4  (9-2) 

where  C„3  is  a  characteristic  constant  of  the  medium.   The 
data  show  that  Eq.  (9.2)  applies  to  air  or  steam  up  to  about 
10  atmospheres  pressure  with 

'      CB3  =  9.49  x  10 3  m2  °K/cmb  (9.3) 

In  this  study  we  assume  that  Eq.  (9.2)  can  be  extrapolated 
up  to  about  100  atmospheres  pressure.   We  also  assume  that  an 
expression  of  the  same  form  applies  to  other  media  besides  air 
or  steam  but  that  each  medium  has  its  own  characteristic  value 
of  the  breakdown  constant  C_ . 

a 

In  the  two-fluid  system,  the  separate  values  of  the 
breakdown  constants   CR,   and  C„2   are  usually  known,  but  the 
breakdown  constant   CR3   of  the  resulting  mixture  is  seldom 
known.   In  the  absence  of  adequate  test  data  bearing  on  this 
point,  we  tentatively  assume  that  CR3   can  be  estimated  from 
the  hypothetical  relation 

CB3  "  k;    [xRlCBl  +  (1  "  X)R2CB2]      (9*4) 

Upon  substituting  Eq.  (9.2)  into  (9.1)  and  dividing 
through  by  the  mass  flow  rate,  we  find  the  electrical  work  out- 
put per  unit  mass  of  fluid  in  the  form: 


2   2 
P     /£  C _-  R-\ 

f  •    (  i         )    H   =  CP3  (Ts4  "  Ts5  >  <9'5> 
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By  use  of  standard  perfect  gas  relations,  Eq.  (9.5)  can 
be  changed  to  a  more  useful  form.   In  this  connection  set 

(Ts4  "  Ts5>  =  ATs  (9'6) 

and  2 

a  CB3  Po  (9.7) 

3 

T 
o 

2 

where  P   =  ambient  pressure,  N/m 

T  =  ambient  temperature,  °K 

Notice  that  3   represents  a  dimensionless  version  of 
the  electrical  breakdown  constant. 

Eq.  (9.5)  now  yields  the  important  result 


^S    =     6     fl3li\f!5ilf!o    ,  ■:  ^_-_    ;,, 


Ts4  1\P^    \T 

The  corresponding  gross  electrical  work  output  per  unit 
mass  of  fluid  is  then  simply 

W*  =  C„0  AT  (9.9) 

e     P3    s 

Unfortunately,  the  known  value  of   e   and  the  typical 
experimentally  measured  values  of   C~   are  very  small.   Con- 
sequently, the  dimensionless  breakdown  constant   3   is  ex- 
tremely small  compared  to  unity.   Hence  the  drop  in  stagnation 
temperature   AT   and  the  gross  electrical  work  per  unit  mass 
W*   are  also  typically  very  small.   This  is  a  very  basic  and 
serious  limitation  on  the  performance  that  can  be  obtained  from 
an  EHD  generator. 

Treating  the  process  4+5  as  isentropic,  we  may  write 
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Y3 

(y3  -  i) 

P        =   P        [    S5)  (9.10) 

By   analogy  with  Eq.     (4.8)    we  write   the   continuity   rela- 
tion  as 

K   =  fMA)    =  S5      5    -   f(M,.)    =   m_  (9.1D 

/*3   Ts4^3  /*3   Ts5/Y3 

It  then   follows    that 

f(M=' =7Sfe)f<v 

This    fixes      f(M_)      whereupon      M_      may  be    found  in    the 
usual  way    from  Eqs .     (4.16),    (4.17)    and    (4.18).      However,    con- 
vergence  should  now  be   very   rapid  since  we   may    take    as    a 
first   approximation 

M_  w   M.  (9.13) 

5  4 

This  solution  now  suffices  to  fix  all  properties  at 
states   4   and  s4  . 


41 


10.      Condenser/Separator 

The   primary    fluid  which    leaves    the    condenser/separator 

at  station   6    is    assumed   to  be    compressed   liquid   at  known 

static  pressure      Pfi    =   P_    =   P.       and   at  known   static    temperature 

T.    =  T_    .      The   kinetic  energy  of   the    liquid   at   station   6    is 
b  A 

assumed   to  be   negligible. 

Temperature      Tfi      must  be   equal    to   or   greater   than    the 
ambient   temperature      T        in   order   to   satisfy    the    requirements 
of   a  heat  balance   on    the    condenser/separator.      No   attempt  is 
made    in    this    analysis    to    formulate    this    actual  heat  balance; 
instead     Tfi      is    simply    treated   as    a   given   or  known  quantity. 
However,    it  may  be   of   interest   later   to   investigate    the 
effects   of   varying     Tfi      on    the   overall  performance   of   the   EHD 
system. 

The   secondary    fluid  which   leaves    the    condenser/separator 
and  enters    the   ejector  at  station    2    consists   primarily   of   a 
low  molecular  weight  noncondensing   gas.      It  also    contains    some 
secondary   vapor ,   but   the  mass    fraction   of    this    vapor   is    treated 
as   negligible    in    this    analysis.      Thus    the    secondary    fluid   at 
station   2    is    treated  as    dry   gas.      The   kinetic   energy   at  station 
2    is   not  negligible    and  must  be    taken   into    account  in    the   ana- 
lysis . 

In   analyzing   the    thermodynamic  process    in    the    condenser/ 
separator,    it   is    again   useful   to   employ    the   concept  of  avail- 
ability  in   steady    flow.      The   availability    function   is   always 
defined  with   respect   to    some    suitable    reference    state.      Nor- 
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mally  this  is  fixed  by  the  conditions   P   ,  T   of  the  ambient 

o  o 

atmosphere.      In   some    cases,    however,   we   seek   to   characterize 

the   performance   of   some   component  of   the    system  in   a  manner 

which   is    independent  of    the   ambient  conditions      P      .    T      .      In 
r  o  o 

such   cases    it  is    useful    to   choose    as    a   reference   state   some 
condition    that   actually  occurs    in    the    component  in  question. 
We   used  such   a  procedure    to    good   advantage   in    analyzing   the 
diffuser   and  in   defining   the    diffuser  effectiveness      nD    .      We 
now   use   it  again   in   analyzing   the   condenser/separator   and  in 
defining  a    friction    loss    coefficient      c-     which   characterizes 
the   operation  of   this   component.      For   this   purpose  we    choose 
as    the    appropriate    reference   state   of   zero   availability    the 
static   conditions      Pg    ,    T,      in   the  condenser/separator. 

It  also    clarifies   matters    further   if  we    choose    the    con- 
dition     P_    ,    Tr      as    the   datum  state   of    zero  enthalpy   and    zero 
6  o 

entropy.      In   this  way  we   ensure    that  our  analysis    in    this   sec- 
tion makes   no   reference  whatever   to   the   state      P      ,    T 

o  o 

Remember  also    that     P-   =  P-   =   P,      and      T,    =   T_    . 

6^1  62 

Thus  we   may  now  write 

s5  =  sm+  cp3  zn(^j  -  r3  ln(^j  (10.1) 

S2    =    CP2    *n(*f)    "    R2^)   =   °  il0'2) 

s„    =    0  (10.3) 

6 

The  quantity  s  in  Eq.  (10.1)  is  the  same  as  previously 
defined  in  Eq.  (5.5).  It  expresses  the  entropy  increase  asso- 
ciated with   mass   mixing. 
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The   availabilities    at  stations    5,    2   and  6    are    therefore 


ip-    =    C    ,     (T    _    -    T0)     -    Tc    sc 
r5  p3        s5  2  6      5 


*2    =   Cp2     (Ts5    -    V    "    ° 


*6    =  h6    "    T6    S6    =   ° 


(10.4) 
(10.5) 
(10.6) 


The   irreversibility   I    of   the   process  which   occurs    in    the 
condenser/separator  with    respect   to    the   specified  reference 
state      Pfi    ,    Tfi      may  be   expressed   in    terms   of    the   above    avail- 


abilities   as    follows 


I    -    1|>-    -     (1    -    X)     *,    -    x>^ 


(10.7) 


We  postulate  that  this  irreversibility  can  also  be  corre- 
lated with  the  inlet  kinetic  energy  in  terms  of  an  empirical 
friction  loss  coefficient   c-   as  follows 


V 


I  =  c 


where    as   usual 


/Cf   ^3   R3    Ts5 


_^5 

(Ts5/T5> 


(10.8) 


s5 


1   + 


<Y3   -    1> 


Mr 


(10.9) 


Of   course    the   above    friction   loss    coefficient      cf     must 
be   estimated   from  suitable   experimental   data. 

Assuming   that      c-      is    a   specified   constant,    the   trial 
solution  obtained   for   an    arbitrary    trial   value   of     x      and 
arbitrary    fixed  values   of   the   various   other  input  parameters 
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will  not  necessarily  satisfy  Eq.  (10.7).   It  is  therefore 
convenient  to  represent  the  residual  error  in  this  equation 
in  the  form  below.   The  factor   R-,  Tc      is  inserted  on  the  left 
to  make  E(x)   dimensionless .   Thus 

R3  T6  E(x)  =  4»5  -  (1  -  x)  m2   -  0  -  I         (10.10) 

The  foregoing  relations  can  be  assembled  and  rearranged 
as  follows.   Let 

Y3 


(Y,  "  1) 


(10.11) 


K^sfUf-1"  (10-12) 

Then    the    final   energy   balance   equation   reduces    to 

E(x)    =   K      +   K,    x  (10.13) 

o  1 

By  holding  all  other  input  parameters  fixed  and  allowing 

x   to  take  on  a  range  of  possible  values,  it  will  normally  be 

possible  to  find  by  trial  a  value  of   x   such  that 

E(x)  =  0  (10.14) 

This  then  is  the  value  of   x  which  is  consistent  with 
the  other  specified  input  parameters  and  which  defines  a  pos- 
sible EHD  thermodynamic  cycle. 
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Of  course,  if  the  various  fixed  input  parameters  are  not 
properly  chosen,  Eq.  (10.14)  cannot  be  satisfied  for  any  trial 
value  of  x  which  lies  in  the  useful  range   0<x<l   and  no 
solution  exists  under  these  conditions. 

If  a  solution  exists  it  can  be  found  by  some  suitable 

iteration  technique.   One  useful  method  is  to  set  E(x)  =  0 

in  Eq.  (10.13)  and  then  solve  this  equation  for   x  but  to 

relabel  this  as   x     .   Then 

new 

Assuming  that  conditions  have  been  so  specified  that  a 
solution  for   x   actually  exists,  and  assuming  that  an  initial 
trial  value  of   x  has  been  chosen  that  is  sufficiently  close 
to  the  true  root,  repeated  application  of  Eq.  (10.15)  should 
finally  yield  convergence  to  the  true  root  itself. 

Incidentally,  any  trial  solution  for  which  E(x)<0   re- 
presents a  physical  impossibility  because  it  violates  the 
second  law  of  thermodynamics.   Any  trial  solution  for  which 
E(x)>0   satisfies  the  second  law  but  represents  a  condition 
in  which  the  losses  in  the  condenser/separator  exceed  the 
value  prescribed  by  Eq .  (10.8).   Only  the  case  for  which 
E(x)  =  0   satisfies  both  Eq.  (10.8)  and  the  second  law. 
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11.      Cycle  Efficiency 

Once  a  solution  has  been  found  such  that  E(x)  =  0  ,  it 
becomes  a  straightfoward  matter  to  calculate  the  correspond- 
ing overall   cycle   efficiency. 

Because   of    the   pressure   drop    through    the  boiler/super- 
heater   ,    the  pressure      P_  =  P    -.      at  boiler  inlet  is    slightly 
higher   than    the    stagnation  pressure      P    ,      of   the   primary 

fluid  at    the   ejector   inlet.      We   assume    that   the   ratio      P    .  /P    - 
J  si      s7 

is    a   specified   constant.      Hence  with      P    ,      specified,      P    _, 

c  si   ^  s7 

is  also  known.   Also  the  pressures   P,  =  Pn   are  known. 

D  1 

* 

Consider   the  ideal  pump  work     w        done   per   unit  mass   of 
primary    fluid.      For   the   present   consider   the  hypothetical   case 
of   a   reversible   pump.      Also  note   that  density   changes    across 
the   pump   are   negligible.      Hence  we   may  write 

*         (Ps7   "    V 

w     =  — — —  (11.1) 

P  P6 

The  ideal  gross  electrical  work  output  per  unit  mass  of 

mixed  fluid  has  earlier  been  established  in  Eq.  (9.5).   It  is 

* 

denoted  by  symbol  w 

Let  us  now  denote  the  net  useful  electrical  work  output 

per  unit  mass  of  primary  fluid  by  symbol  w    .   It  may  be 

*  * 

related  to     w        and     w        in   the   following  way 

e  p  '        J 

/  *  *\ 

T\        W  W      \ 

w      .    =     — — - £  Joule/kg  (11.2) 

net        y      x  rip/ 

where      n      ,    the   excitation   efficiency,    allows    for  the   small 
electrical   power  expended   to   excite    the  system  and  where      n       , 
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the  pump  efficiency,  allows  for  irreversibilities  in  the  pump. 

The  heat  input  in  the  boiler/superheater  per  unit  mass 
of  primary  fluid  may  be  written 

* 

«i*-  (hsi  -  v  -!j  (ii-3) 

where      hfi      is    the  enthalpy  of   the   primary   liquid   at   conditions 

P..    ,    T        as    listed  in   suitable    tables   of   properties, 
o  o 

The   overall    thermal   efficiency   of    the    cycle    can   now  be 
calculated    from   the   simple    formula 

nc=^  (11.4, 

Both      n        and     w  are    useful   parameters  which    charac- 

c  net  e 

terize    the  overall    thermodynamic   performance   of   the  EHD  sys- 
tem. 

Eq.     (11.2)    shows    clearly   how   the   ejector   serves    as    a 
kind  of   amplifier  which    increases    the   electrical  work  output 
per  unit  mass    of  primary    fluid.      This   effect  is    shown  by    the 
presence   of   the   parameter     x     in    the    denominator  of    the    first 
term.      Notice    the   beneficial  effect  of   a   low  value   of      x     on 

parameters      w  and     n      .      Unfortunately,    the   value  of      x 

e  net        c  2 

cannot  be  stipulated  independently  in  advance;  it  is  fixed 
by  the  other  specified  input  parameters  as  explained  in  the 

previous  section. 

* 

Notice    that   owing   to    the    fact   that     w        is   normally 

3  e 

very  small,  it  is  possible  in  some  circumstances  for  w   , 

*  '  r  net 
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to    take   on   values  which   are    actually   negative.      When    this 
happens    it  means    that  electrical  power  output  is    insufficient 
to   drive    the   pump. 

The    relationships    developed  in    this    and   the   preceding 
sections   make    it  possible    to   carry   out  systematic  parametric 
studies    of   various   one-fluid   and   two-fluid  EHD   cycles,    with 
realistic  allowances    for   the   various    losses    that  occur.      Such 
studies    can   establish  optimum  design   parameters    and  perfor- 
mance   limits    under  various    circumstances.      These    results    in 
turn   can    finally  permit  informed   conclusions    to  be   drawn   con- 
cerning  the   ultimate    feasbility   of    this    general    type   of  EHD 
power  generator. 
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12.   Calculation  Sequence  for  Two-Fluid  System 

In  this  section  we  summarize  the  various  equations  pre- 
viously derived  in  the  approximate  order  in  which  they  would 
be  used  in  the  calculation  of  system  performance.   Also  listed 
are  the  initial  input  parameters  whose  values  must  be  specified 
in  order  to  start  the  calculation  and  various  further  input 
parameters  whose  values  must  be  specified  in  order  to  proceed 
with  various  successive  stages  of  the  calculation. 

Initial  Input  Data:   px  =  P2  =  P6  '  Tl  '  Psl  '  Tsl  '  (hsl  "  hl}  ' 

R  =  8315   Joule/kg°K 

*i  "   I1  -IHTp^7pp-)  (12-1) 


<hsl  -  hl> 


Ri  =  (^T-|cp1  U2-3) 


"i"^  <12-4) 


Further   Input   Data:      W      ,    W_   =  W~    ,    T0    =   Tc    ,    P    (T^) 
c v    '      G  2    '      2  6    '      v      6 


"    -      %       IP2    -    PV(T6) ]  <12-5> 

Verify    that      oo<<1 
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V, 


Further  Input  Data:   y2  /  I  77- )  =  Y 


Ml  =  /(Yl  -  1) 


si 


-  1 


(12.6) 


M, 


M  Ai  Ri  Tsi 


M. 


l/V  Y2    *2    x2        V/Tsl/T1 


(12.7) 


Verify    that      M2<1 


(12.8) 


¥)  -  -  P^  »22 


Ts2    =    T2(Ts2/T2} 


(12.9) 
(12.10) 


y2 


s2 


s2 


(12..  11) 


Ps2    "   P2      P 


s2 


(12.12) 


Further  Input  Data:      x      (trial   value) 


R2   "  W. 


R3   =   xRx   +    (1    -    x)     R2 


(12.13) 
(12.14) 


y2    R2 

p2  (Y,    -    1) 


(12.15) 


C    ,   =   X   C   ,    +    (1-x)     C0 
p3  pi  p2 


(12.16) 
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*p3 


"3         (Cp3    "    R3) 


(12.17) 


?    o    =   tt1-    [x   C    ,    T    .    +    (1    -    x)     C    0    T    0] 
s3        C   3  pi      sl  p2      s2J 


(12.18) 


f(M1) 


=    M, 


Yl    "    1  2 

1    +      ^-      Ml2 


(Yjl   +   1) 

2(yx  -  i: 


(12.19) 


f(M2)     =    M2 


Y2    "    1  2 

1    +      -i-2—      M22 


(Y2    +   1) 
2^2    "    1) 


(12.20) 


Ps2\     /Ri    Tsl   ^2    ^(M2> 


A2  (1   "    X)    \Psl/\/R2    Ts2    Yl    f(Ml} 


(12.21) 


(Ax/A2) 


A3     (1  +  A1/A2) 


(12.22) 


Further  Input  Data:   P   ,  T 
c o    o 


si  -  cPi  £n  (W  -  Ri  ^(r1 


(12.23) 


s2    "   CP2    *n   (ff )  "    R2    £n  (Jf 


(12.24) 


to(5?)-5rjcp3  ^iM"  xsi~  (1"x)  s2 


(12.25) 


£n    (P      /P    ) 
sx'    o 


P        =   P      e 
sx  o 


(12.26) 
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(1    -    A  /A3)     /P    2\       R3    T    3    yV 


(Y3  +    1) 


'max        \y-5   +    1 


2(Y,    "    1) 


(3    -    y,) 


F1  (M   )     =    f 
n' 


^'  4'^  «.f VJ " l 
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(12.2  8) 


Verify    that      f(M  )    =    f<f  (12.29) 

2  x  max 

For   first  approximation   set     M     =  M~  (12.30) 

Iterate   using  Eqs.     (12.31),    (12.32),     (12.33). 

(Y3  +  i> 


(  (*3   ~    X)  2\  3 

F(M    >     =    f       1    +     X^-^ '     M„  -    M  (12.31) 

n  I  z  n    I  n 


-    1) 

(12.32) 


F(M   ) 

M(n  +   1)    "   Mn   "   F^THT  (12-33) 

n 

Iteration  converges  to  fix  M 

A 


Ml      \       Yl  '  X     2 
g(M  )  =  ± *_  \/l  +  (-^5 )  M/  (12.34) 

(1  +  Yl  M i  )   V         * 


M2      .       Y2  "  1 
g(M  )  =  1 \  i  +  (_fL )  m/  (12.35) 

2     (1  +  Y2  *V  ) 


v-1 


g(My) 


=  g  = 


R3  Ts3 


x 


,R1  Tsl   +  (1  -  X) 


Y3   I  g(Mx)  /  y1  g(M2)   V   y2 


R  T 


s2 


(12.36) 


M 


2  = 


(1  -  2y3  g  )  -  y  1  -  2  (y3  +  1)  g 
1  -  Y3  (1  -  2y3  g2) 


(12.37) 


Eq.  (12.37)  applies  provided  that  g  ^  g^  = 


(Y3  "  1) 


(12.38) 


If   g  =  gOT   the  solution  reduces  to 


MY  = 


(Y3  "  1) 


2Y 


(12.39) 


Once  M   is  known  from  Eq.  (12.37)  or  (12.39! 


l) 

(Y3  +  1) 

f(M  )  =  M 

y    y 

(       fr," 

1+     2 
V 

"J 

f(M  ) 

p   -  p     x 

sy     sx  f(M  ) 

Further  Input 

.  Data:   rir 

? 

(12.40) 


(12.41) 


l/n. 


p  ^  =  p 


sy 


s3    sx  I  p 


sx, 


(12.42) 


f(M3)  =   - 


sx 


s3, 


f(M  ) 
x 


Verify  that  f(M3)  =  f<f 


Max 


(12.43) 
(12.44) 


54 


For   a   first  approximation   take      M_   =  M 
r^  3  y 


(12.45) 


Iterate  using  Eqs .  (12.31);  (12.32),  (12.33) 
Iteration  converges  to  fix  M~  . 

T 


s3 


(Y3  "  X)    2 
1  +  L^_'  M32 


(12.46) 


T3   =    Ts3/(T83/T3) 


(12.47) 


Y 


(Yt  -  i) 


s3 


s3 


(12.48) 


P3   =   Ps3/(Ps3/P3) 


Further  Input   Data:      r\       ,  I  — 


P    a    =   p    ->   e 
s4  s3 


1  -   n 


(12.49) 


(12.50) 


f(M4)     -    f    -l!sl)(^|f(M3) 


(12.51) 


Verify    that      f(M4)    =    f<fmax 


For  a   first  approximation   take     M.    =    [»•— 1    Mo 

Iterate   using  Eqs.    (12.31),     (12.32),     (12.33) 

Iteration   converges    to    fix     M,    . 

4 


(12.52) 
(12.53) 


T       A       ~      T       "5 

s4  s3 


(12.54) 
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Ts4  fY3    "    X)  2 

^  =    1    +    V— =2— g — : '     M/  (12.55) 


T4    =    Ts4/(Ts4/T4)  (12.56) 

Y3 


P  T       (Y3    "    1] 

s4  s4 

p^  =      ^-  (12.57) 

F4  x4 

?4  =  Ps4/(Ps4/P4)  '  <12-58> 

-12 
Further  Input  Data:   £  =  8.854  x  10    farad/m,  C-,  ,  C_.~, 

Dl      QZ 

3   2 
(C_  =  9.5  x  10  m   °K/cmb  for  air  or  steam) 

CB3  "  SJ  [X  Rl  CB1  +  (1  "  X)  R2  CB2]  <12-59' 


£    CB32    Po 


(12.60) 


ATs    =  3  (Irff^iV^2  fx  +  ,!L^,.  2f /<Y3  "aL, 


Ts4  \    2^3      i\Po  ilTs4' 

AT 

S5      =    (1   -   =-£)  (12.62) 

T    .  s4 

s4 


T    -      T3/(Y3    "    1J 

Ps5      "   Ps4    't5?  <12-63' 

s4 
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f  ( Me )    =    f   =       /-^-  I  -^i  1 

5  vTs4  yS5i     4  <12-64> 

Verify    that      £(M„)    =    f<f     _  (12.65) 

For   a   first  approximation   take      M,.    =   M.  (12.66) 

Iterate   using  Eqs.   (12.31),     (12.32),     (12.33) 
Iteration   converges    to    fix     M_    . 

sf  J  ■ 1  +  ir-T-j  Ms2  (12-67) 

T5   ■   WWV  (12.68) 


<y3  -  n 


Ps5\      /Ts5 


P5    /       \T5     /  (12.69) 

Further   Input  Data:       cf 

sm  =    R3    in   R3   -    x    Rx    £.n    (x   1^)  (12.71) 

-    (1   -    x)     R2    in    [(1   -    X)     R2] 

«» = -  (Jf )  *  ^MH  -  ^- -  fr 

Sm  Cf^/T5\M2 

"    R"     "   — 2—  (fj  )M5  (12.72) 

Kl   =    R^(t^  "    y  C12.73) 
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Verify    that      K    <0 

2  o 


X 


K 

new  K. 


(12.74) 
(12.75) 


Iterate    from  Eq.      (12.14)    until      x      converges. 

Further   Input   Data:       (psl/ps7)     /    P6    >    (hgl   -   hg)     ,    nx    ,    n 


Ps7   "   Psl/(Psl/Ps7) 


e  p  J        s 


w 


*   =    (Ps7    "    V 


(12.76) 
(12.77) 

(12.78) 


* 
n      w  w 

w      .    -I-5-S  -  JEL 


W 


q,n    =    (h    ,    -    hj     -    -J 
^m  si  6  n 


RESULT  (12.79) 


(12.80) 


w 


net 


C         q 


in 


RESULT  (12.81) 
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13.   OPTIMUM  VELOCITY  RATIO 

The  calculation  sequence  summarized  in  the  preceding  sec- 
tion reveals  that  a  comple  cycle  calculation  requires  us  to 
specify  numerical  values  for  about  thirty  input  parameters. 
One  of  these  is  the  velocity  ratio  (V../V2)  =  y  .   The  calcula- 
tion procedure  outlined  above  then  permits  us  to  find  the  cor- 
responding value  of  the  mass  flow  ratio   x  ,  assuming  that 
there  exists  a  value  of   x  which  satisfies  the  governing 
thermodynamic  relations  for  the  specified  values  of  the  other 
input  parameters. 

If  we  now  change  the  value  of  y  but  hold  all  other  in- 
put parameters  fixed,  we  can  repeat  the  above  procedure  and 
find  a  corresponding  new  value  of   x  ,  if  such  a  value  exists. 
It  is  evident,  therefore,  that  x  becomes  some  definite 
function  of  y   as  long  as  we  remain  in  the  domain  where  a 
real  solution  exists.   Of  course  all  other  dependent  variables 
of  the  cycle  are  also  functions  of  y  .   In  particular,  the 
overall  cycle  efficiency  n    is  some  definite  function  of  y  . 
Moreover,  for  prescribed  values  of  the  other  input  parameters, 
there  will  be  some  definite  value  of  y  ,  let  us  designate  it 
as  the  optimum  value  y     ,   which  yields  the  greatest  value 
of   ti    that  is  possible  under  the  specified  input  conditions. 
Our  basic  purpose  therefore  is  to  determine  for  any  prescribed 
values  of  the  other  input  parameters,  the  value  of  y   .   an(3  the 

correspondina  value  of   ri  )_   •   Of  course  the  values  of   x,,  . 
c  c  max  opt 

and  other  dependent  variables  are  also  of  interest. 
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Once  a  working  computer  program  is  available  to  establish 

values  of   y   ,  ,   i"i  )  „v  ,      x     and  so  on  for  prescribed 
opt     c  max     cpt 

values  of  the  other  input  parameters,  we  can  then  proceed  to 
study  the  effects  of  changes  in  these  other  parameters  and  to 
search  systematically  for  such  values  of  the  most  important 
input  parameters  as  will  yield  the  best  overall  performance 
of  the  system. 
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14.   CALCULATION  SEQUENCE  FOR  ONE-FLUID  SYSTEM 

The  calculation  sequence  summarized  in  Section  12  simpli- 
fies for  the  special  case  where  the  same  condensible  fluid  is 
used  for  both  the  primary  and  secondary  streams.   This  special 
case  is  summarized  below. 

Initial  Input  Data:   P,  =  P„  =  Vr    ,   z-,  ,   Pi/   T-  =  T£ 
s 1     2     6      1      si      2     6 

Assuming  pressure   P,  *  and  quality   z,   are  specified,  we 
may  find  from  tables  of  properties  of  the  saturated  fluid 

hl  =  hg  "   t1  -   ^fg  (14.1) 


Sl   =    Sg    "    (1    "     Vsfg  (14*2) 

Then  with      s,       fixed   and     P    ,       specified,    we  may   also    find 

T  and     h    ,      from  tables   of  properties    of   the   superheated 

vapor. 

In  addition  with   Tfi   specified  somewhere  in  the  range 
T   <_  Tr  <_  T,  ,  and  with  P^  =  P,   specified,  we  may  find  ¥l-    f    the 
enthalpy  of  the  liquid  at  state  6  from  the  tables  of  properties. 

The  primary  fluid  can  now. be  approximated  by  an  "equivalent 
perfect  gas"  whose  properties  are  as  follows 

-1 
f  2,n(T   /T  )\ 

(hsl  "  V 
CP  =  (Tsl  -  Tp  (14'4> 

R  =  (Y  "  1)C  (14.5) 

Y     P 
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The  secondary  fluid  is  represented  by  the  same  "equivalent 
perfect  gas"  and  is  therefore  assigned  the  same  values  of 


Y  ,   C    and   R 
P 


Further   Input   Data;       (V^/V   )    =   y 


M. 


(Y    -    1) 


'I    ) 


-    1 


(14.6) 


n. 


vi  i  Ti    • 

-2)  k  Mi 


Verify   that     M2    <    1 


(14.7) 
(14.8) 


s2 


=    1   + 


(^) 


M 


s2       =    T2(Ts2/T2) 


(14.9) 


(14.10) 


pf      "    <Ts2/T2> 


Y/(Y    "    1) 


(14.11) 


Ps2      "   P2(Ps2/P2} 


(14.12) 


Further  Input  Data:   x  (trial  value) 


T  -j   =  xT,  +  (l-x)T0 
s3       si  s2 


(14.13) 


f(Mx)  =  M± 


f  (M2)  =  M2 


i  +  (:l-t^)m21 


-(Y  +  1)/2(y  "  1) 


_-(Y  +  D/2(Y  "  1) 


1  +  (Y  2  1)M2 


(14.14) 


(14.15: 


Al^  _fPs2\JTsl    f(M2} 


x 


PsljVTs2    f<Ml>        tl   -   x) 


(14.16) 
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a 


Ax/A2 


(1    +   A1/A2) 


(14.17) 


Further   Input   Data:      P       ,      T         (ambient   conditions) 
£ o  o 


=   cp£n 


o  v      o/ 


(14.18) 


88  Vn   r 


-    R£n 


(14.19) 


£n 


sx 


R    IS    "n 


's3 


x   s,    -    (1   -   x)     s2 


sx 


P      e 
o 


£n(P      /P    ) 
sx'    o; 


(14.20) 


(14.21) 


f(M   ) 

x 


max 


(1   -   AVA,)  /P    0\  .    T    . 

32\\l_sl     f(M2) 


(1  -   x) 


sx/    «     s2 

1)/2(y    "    1) 


\y  +  i 


Verify   that      f(M  )    =    f   <    f 

X  iuci..X 


For   a   first  approximation   set     M      ~   IA0      and   iterate 
Iterate   using   Sqs.     (14.26),     (14.27),     (14.28). 


(14.22) 

(14.23) 

(14.24) 
(14.25) 


F(Mn)  =    f       1    +     ( 


Y 


1      2\(y   +    D/2(Y    "    1) 


-  -)  M' 
2  n 


-   M 


n 


(14.26) 


F' (M    )       =    f 
n 


<^4  ♦  iH±<)13 ' Y>/2(Y "  "i   <14-27) 


L_    m 


=   M      - 


F(M    ) 
n 


(n   +    1)  n        F'  (Mn) 


(14.28) 
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Iteration  converges  to  fix  M 


M. 


=   ^-       s/l  +    (I^-ijM^ 

(1  +  yMJ)  ^  x 


(14.29) 


M, 


g(M2) 


(1   +    YMJ) 


A   +    (I^-i)M^ 


g(M-  )       =    g    =    I 

2 


x     "\/Tsl    ,    (1  -  xn/Ts2 
g(M1)   Vts3        g(M2)        ill   3 


M 


_  (1  -  2yg2)  -   y/l  -"  2(Y  +  1)  g2 


(14.30) 


Y  1  -  Yd  "  2Yg  ) 

It  can  also  be  shown  that  in  the  special  case  where 


(14.31) 


(14.32) 


1   /v  -  l 
g  =  g^  =  —  v-i — 2^ /  tne  above  solution  degenerates  to 


Once  M   is  found  from  Eq.  (14.32)  or  (14.33) 

y 


f  (My)   =  My 


sy 


=  P 


1+  (I^)My 
f  (Mx) 


-  (y  +  1)/2(Y  -  1) 


sx  f(M  ) 

y 


Further  Input  Data:   n.., 


(14.33) 

(14.34) 
(14.35) 


1/tU 


s3 


=  P 


sy 


SX   P 


sx 


f(M3)  =fe) f(MX) 


(14.36) 


(14.37) 
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Verify  that   f(M-)  ■  £  <  f  (14.38) 

o         max 

For  a  first  approximation,  take   I'  Z   M  (14.39) 

•5   y 

Iterate  using  Eqs .  (14.26)  ,  (14.27)  ,  (14.28)  . 
Iteration  converges  to  fix  1^  . 


Ts3  =    .   (y  -  1)  2  •  (14.40) 

T  2    u3 

T,   =  T  ~/-(T  ~/T,)  (14.41) 


pq,         y/(y  -  l) 

gSLi  =  (1'   /TO  (14.42) 

P3   =  Ps3/(Ps3/P3)  (14'43) 

A- 


l3 

Furtner    Input   Data:       ru    /       * — 
u  a4 

-(1   -    nn)  (Y    M?/2) 
Ps4      -   Ps3    e  (14.44) 

fU4)=f    -    (fS)(^)    f(M3)  (14*45) 

Verify  that  f(M,)  =  f  <  f  (14.46) 

1  4'         max 

For  a  first  approximation  take   M,  Z   |  ^-j  M-,  .       (14.47) 

Iterate  using  Eqs.   (14.26),  (14.27),  (14.28)  . 

Iteration  coverges  to  fix  M..     . 

T  .   =  T  .  (14.48) 

S4      s  J 


T 


;4 


1  +  (Y  -  1)  m2 
T ,       X  +     2    ^4 
4 


(14.49) 
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T4      =  Ts4'/(TS4/T4)  (14.50) 


^          =    (T34/T4)Y/(Y    1}  (14.51) 

P4      "  Ps4/(Ps4/P4>  (14-52) 

Note:   If  it  be  desired  to  omit  the  diffuser,  merely  set 

nn  =  1   and   At/A.  =  1  .   The  result  is 
D  o   4 

Ts4     =  Ts3  (14.53) 

T4      =  T3  (14.54) 

Ps4     =  Ps3  (14.55) 

P4      =  P3  (14.56) 

M4      =  M3   etc  (14.57) 

and  all  effects  of  a  diffuser  disappear  from  the  result. 

-12 

Further  Input  Data:   e  =  8.3b4  x  10  farad/m   for  all  gases 

Cg=  9 . 5  x. 10   m   K/cmb   for  air  or  steam 


2 
ECBP0 

T  2 
o 


(14.58) 


AT 
s 

s4 


^  (y  fe)2 11  +  ^  41/<Y^4-59' 


T  -     =  T  ,  -  AT  (14.60) 

so        s4      s 

Psi     =  Ps4(Ts5/Ts4)  (14.61) 


f('v    =(fe)(fe?)f(-v  (14'62) 
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Verify  that   f(MK)  =  f  <  f  (14  g31 

3         max  UH,DJ^ 

For  a  first  approximation  take   M-  ~  M,  .  (14.64) 

Iterate  using  Eqs.  (14.26),  (14.27),  (14.28)  . 
Iteration  converges  to  fix  M,-  . 


3rt 


T,- 


1  +   (y  ~    $     M2 

1        2      o  (14.65) 

=      Ts5/(Tso/T5)  (14.66) 


(5f)         -    (Ts5/T5)Y/(Y    "    "  <14'67> 


(14.68) 


o  sr        sb/     5 

Further  Input  Data:   c- 

Note:   For  a  single  fluid  cycle   s_  =  0 

3  ■*      m 

Ki    =TT^TT  (ff "  x)  (14-70) 

E(x)    =  KQ  +  K,  x  (14.71) 

K 

o 

Xnew     "  K,  (14.72) 

Return  to  Eg.  (14.13)  an(j  iterate  until   x   converges  and 
E(x)  =0  . 
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Further  Input  Data:   (psl/p37)  t      P6  /   (hsl  "  hg)  ,   n   ,   n 


s7 


"  ^l^sl^s?* 


w 


=  C   AT 

P    s 


P 
(14.73) 

(14.74) 


w 


(Ps7  "  V 


(14.75) 


w 


net 


n  w     w* 
x  e      p 

x    n7 


RESULT  (Joule/kg)  (14.76) 


lin 


si    6'    n 

P 


(14.77) 


w 


net 


qin 


RESULT 


(14.78)- 
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15.   SAMPLE  CALCULATION  FOR  A  ONE  -^LUID  SYSTEM 

In  this  section  we  illustrate  numerically  how  the  calcu- 
lations for  a  one -fluid  system  may  be  carried  out  for  one 
arbitrary  but  typical  operating  point.   To  find  an  optimum 
operating  condition  would  require  the  repetition  of  similar 
calculations  for  a  large  number  of  such  trial  operating  points, 
No  attempt  is  made  in  this  section  to  find  such  an  optimum.   * 

We  arbitrarily  choose  to  calculate  a  case  which  is  in 
some  respects  similar  to  that  identified  as  case  B  in  the 
appendix  of  Ref.  (1) .   Thus,  we  have  for  both  of  these  cases  - 
Medium 
Maximum  pressure 

Quality  at  state  1 

Ejector  effectiveness 

Friction  factor 

Dielectric  strength 

Diffuser  characteristics 
Excitation  efficiency 

Pump  efficiency 

Pressure  ratio  across 
bo  i  le  r 

Nevertheless,  despite  these  similarities,  these  two  cases 
are  not  identical.   Case  B  treats  the  densities  at  states  1,  2, 
x,  y,  3,  4,  5  as  equal.   The  present  analysis  takes  the  true 


Steam 

Psl 

1500  psia 

zl 

0.95 

nE 

0.9 

Cf 

0.1 

CB 

9,4  90   m2  ° 

No  diffuser 

nx 

1 

np 

0.9 

Psl/ps7 

1 
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density  changes  at  these  states  fully  into  account.   Also  it 

is  expedient  in  the  present  compressible  analysis  to  define 

the  parameters  r\        and   cf   somewhat  differently  than  they 

are  defined  in  Ref  (1) . 

Ref.(l)  also  shows  that  a  comparatively  good  performance  is 

obtained  by  choosing 

P,  =  P.,  =  Va   =  600  psia  (15.1) 

12     0^ 

and 

tt  =   y  =  4'25  (15.2) 

v2 

Therefore  we  arbitrarily  choose  the  same  values  for  the 

present  numerical  example,  despite  the  fact  that  it  is  not 

entirely  equivalent  to  case  B  of  Ref  (1) .   We  also  choose  to 

set  T6  =  T]_  (15.3) 

The    above    conditions   now   suffice    to    fix   the   following 

properties.      These  are   obtained   from  Keenan    and  Kayes    steam 

tables    and   are    listed  on   page    10A  of   Ref.     (1)  . 

T    ,  =      1115.9    °R 

si 

T,  =      T2    =    T6    =    946.2    °R  (15.3) 


(h    ,    -   h1)    =      78.9   Btu/lbm 
(h    ,    -  hg)    =      673.9   Btu/lbm 

We    designate    ambient   conditions    as    follows 

P      =      14.696    psia   =    1.013   x    105   N/m2  (15.4) 

o  c 

T      =      520    °R  =    2  38.9    °K 

o 
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The  properties  of  the  "equivalent  perfect  gas"  may  now  be 
computed  as  follows 


&n«P8l/Pl) 


(l    *n(1115. 9/946. 2)A     ,  ,10-  ,__  -. 

"   \X £n(1500/600)   j   "  1*2196  (15*5) 

c  m    (hsl  "  V  _       78.9 


p    (Tgl  -  Tx)    (1115.9  -  946.2) 

=  0.4649  Btu/lbm  (15.6) 

r  =  (Y  -  1)  C   -  °-2196  (0  4649) 
K       Y      P    1.2196  ^u-4b^ 

=  0.08371  Btu/lbm  °R  =  65.15  ft  lbf/lbm  °R         (15.7) 
The  following  constants  will  prove  useful  for  later  appli- 
cation 

Y  "  1   =    0.1098  , — ^-r,    =  4.5537 

2  (y  -  1) 

Hr1  =   1-1098  iY±l)=    5.0537  (15.8) 

2  (y-1) 

-* =  5.5537  (3  -  y) 


(Y-D   --—  2(Y-i)  =  4'0537 

Also  the  dielectric  breakdown  constant  can  be  expressed  in 
dimensionless  form  as 

2  2 

Q         ZC   BPo    (8.854  x  10"12)  (9.49  x  103)   (1.013  x  10+5) 

p 2 5 

T  (288. 9)Z 

o 

=  0.9678(10)"3  (15.9) 
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The  Mach   numbers   at   the   ejector    inlet   are 


Ml   = 


(Y    ~    1)        T 


(aSi-X)    = 


0.2196     x94.62 


1115.9 


-    1) 


=      1.2781 


(15.10) 


M0    =       ( 


V2        \ 


M,     = 


1         4.25 


(1)  (1.2781) 


=    0.3007 


(15.11) 


The    following    functions   are   calculated   for   later   use 

-    (Y    +    1)/2(Y    -    1) 
f  (M1)    =   Mx(l+     T    2    XM1Z) 

=12781(1   +    0. 1098(1. 2781)2)"    5-50537 


=    0.5552 


(15.12) 


f(M2)    =   M    (14     I-Z-lM22)-(y    +    1)/2(y    -'1: 


=    0.3007(1  +»    0. 1098(0. 3007)2)"5*0537 


=    0.2861 


(15.13) 


max 


■  (rr-r) 

=    0.5907 


2  (Y    +    D/2(Y    "    1) 


5.0537 


=    ( 


2.2196 


g(Mx)   = 


M     ^  1    +    Y    -    1    M    2 


(1    +    Y    M^) 


1.2781     VI    +    0. 1098(1. 2781)2 
(1    +    1. 2196(1. 2781)2) 


(15.14) 


=    0.4639 


(15.15! 


g(M2)    = 


(l  +  ym22) 


72 


=   0.3007  \/   1  +  0.1098  (0.3Q07)2 
(1  +  1. 2196(0. 3007)2) 

=   0.2722  (15.16) 

Stagnation  conditions  at  station  2  are  calculated  next. 


(Ts2/T1)  =  1  +  (Y  ~  1)M22  =  1  +  0. 1098(0. 3007)2 

=  1.0099  (15.17) 

T2  =  T2(Ts2/T2)  =  94^. 2(1. 0099) 

=  955.6  °R  (15.18) 

(Ps2/P2}  =  (Ts2/T2)Y/(Y  "  1}  =  d.0099)5*5537 

=  1.0562  (15.19) 

Ps2  =  Ps(Ps2/P2)  =  600-°  d-0562) 

=  633.7  psia  (15.20) 

In  this  particular  case  the  entropies  at  stations  1  and  2 
happen  to  be  equal.   They  can  be  calculated  either  from  the  stag- 
nation temperatures  and  pressures  or  from  the  static  temperatures 
and  pressures.   The  latter  choice  gives: 

s±   =   s2  =  cp£n(T1/To)  -  R  £n(P1/PQ) 

=   0.4649  £n(946. 2/520.0)  -  0.08371  £n (600 . 0/14 . 696) 
=-0.0322  Btu/lbm  °R  (15.21) 

Further  Data  Input:   x  =  0.472  (trial  value) 

The  above  value  of  x  was  found  by  previous  trial  calcula- 
tions not  shown  here.   The  calculations  which  follow  verify  that 
this  value  of  x  does  indeed  satisfy  the  energy  balance  equation 
of  the  system  to  within  the  round-off  error. 
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Ts3  =  x  Tsl  +  (1  "  x)Ts2  =  °-472(1115-9)  +  0.528(955.6) 
=  1031.3  °R 


x 


(1  -  x)  XP 


p  o 

(A 


T     f(M2) 


sl 


ts2  Hm^T 


(15.21) 


0.472  ,633.7  > 
0.528  l1500.0; 


^F^^^tHH^  =  °'2103    <15-22> 


Ax/A2 


0.2103 
A3    (1  +  A1/A2)  "  1.2103 


0.1738 


(15.23) 


Since   s,  =  s2  in  this  case,  Eq.  (14.20)  simplifies  to 


P       ,        T  - 
o        F     o 


P   =  P 
sx    oe 


0.08371 
4.1876 


[0.4649  g>n(i??1I3)  +  0.0322]  =  4.1876  (15.24) 


=  14.696e 


520.0 
4.1876 


=  967.9  psia 


(15.25] 


(1  -  A,/A-)   P  - 

f<Mx>  "    (1  -  x)    <P^ 

sx 


^if(M2) 
s2 


(1  -  0.1726)  .633. 7> 
0.528      l967.9; 


1031.  3 
955.6 


(0.2861)  =  0.3045  (15.26 


Solving  for   M   from  Eqs.  (4.16),  (4.17),  (4.18)  gives 

2% 


M   =  0.3225 
x 


(15.27) 


For  an  ideal  ejector 


-1 


x 


g(V     U(MX)  |t^3  +  g(Mx)X 


s3 


0.472 
0.4639 


1115.9    0.528 


1031.3         0.2722      U031.3 


955.6 


=    0.3418     (15.2 


M 


(1    -    2y 


2^        \f 
g  )    -  If 


1    -    2(y    +    l)g' 


1   -  Yd   -   2Tg^) 
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Ill    -    4. 


=       [1    -    2.4392(0.3418)^]     -    Vl    -    4  .  4392  (0  .  3418)  _ 
1    -    1.2196    [1    -    2.4392(0.3418) 2] 

=      0.1658  M      =    0.4071  (15.29) 

M 
f(M    )    =      X 


Y  Y    _    i         2    (Y    +   D/2(y    -    1) 

2  y 

=   0^4071 =    0.3716  (15.30) 

[1    +    0.1098(0.1658) ]°-UDJ/ 

f(Mx)  0    3045 

Psy      =   Psx    fWT  =    967--9(8f37lf)    =    793'1   P3ia  (15'31) 


Psv    1/nE  791    1    1/0'9 

Ps3     =  Psx(P£7)  "   967'9    (f6TTT)  =   775*7   Psia  (15*32) 

So 

p 

f(M3)  =  (^)f(Mx)  =  (9^7;9)  (0.3045)  =  0.3799  (15.33) 

s3 

Solving  for  M,  from  Eqs.  (4.16),  (4.17),  (4.18)  gives 

M.   =  0.4183  (15.34) 

In  the  absence  of  a  diffuser  we  may  write 

T  ,   =   T  ,   =  1031.3  °R 
S4       S3 


P  ,   =  P  -,   =  775.7  psia 
s4       s3  c 


M4    =   M3    =  0.4183 


f(M4)   =   f (M3)   =  0.3799 


(15.35) 


The  drop  in  stagnation  temperature  through  the  working 
section  is  given  by  the  following  expressions. 
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AT 
*s4 


,   P  .   T    2 
^Y    ^o    s4 


[1  +  (^-=-^-)M42] 


1/(Y  -  1) 


n  Qg-7o/in\-3/0.1098w775.7  w  520.0. 
0.9678(10)   (TT2T96)  (T4T696)  fmiTT* 


(15.36) 
-3 


[1  +  0. 1098(0. 4183)z] 


4.5537 


=  1.0722(10) 


AT. 


-3 


(15.37) 


=  Ts4(ATs/Ts4)  =  1031.3(1.0722)  (10)  J  =  1.1057  °R 


Notice  that  this  temperature  drop  is  only  about  1.1  °R  or 
about  one  tenth  of  one  percent  of  the  absolute  stagnation  tempera- 
ture! 

Ts5   =  Ts4  "  ATs  =  1031-3  ~  1-1  *  1030. .2  °R  (15.38) 


Ps5   =  Ps4(Ts5/Ts4) 


Y/(Y  -  1) 


=  775  7  (1030-2) 
//D,/  *1031. 3; 


5.5537 


=  771.1  psia 


(15.39) 


f(M5)  = 


Ts5   Ps4 
s4    s5 


1030.2  ,775.7.  (Q    „qq. 

1031.3  (T7T7T) (°'3799) 


=  0.3820 


(15.40) 


Solving  for  M5  from  Eqs.  (4.16),  (4.17),  (4.18)  gives 


M5   =   0.4210 


(15.41) 


(T  _/Tc)  =   1  +  (Y  ~    1)MC2  =  1  +  0. 1098(0. 4210)2  =  1.0195  (15.42) 

S3         D  Z  D 


T5      =      W^.s'V    =  OjM  =   1010'5    °R 


(15.43) 
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To  verify  the  energy  balance  in  the  condenser/separator  we 
calculate 

Ps5  y  Ts5  C*Y    T5 

Ko  =  *n<pf )  -  ttHt  Zn{^  -  -i-(Tf)M5 

„„  ,711.1,         -    cc,-    0     ,1030.2,         0.12196, 1010. 5wn    ,,..,2 
=    £n(600T0)     "    5*5537    An(946.2    } 2 (946.2    }  (°-421Q> 

=   -    0.2330  (15.44) 

T 
Kl  =  (Y  -  1)  (T^"  "  1]    =    5-5537(946?22  "  1]    =  °'4930      (15.45) 

The  final  energy  balance  now  gives 

E(x)  =  KQ  +  K1  x  =  -  0.2330  +  0.4930  (0.472) 

=  -  0.0003  (15.46) 

This  verifies  that  the  energy  balance  is  satisfied  to  within 
round-off  error.   It  also  confirms  that  the  value  x  =  0.472   is 
the  mass  flow  ratio  that  corresponds  to  the  various  independent 
input  conditions  stipulated  in  this  sample  calculation. 

Now  that  the  value  of   x   has  been  confirmed,  we  may  complete 
the  thermodynamic  cycle. 

Since  the  pressure  drop  across  the  boiler  is  taken  as  negli- 
gible in  this  example 

P  ,  =  P  ,  =  1500.0  psia  (15.47) 

s7     si 

The  gross  electrical  work  output  per  unit  mass  is 

w  *  =  C  AT   =  0.4649(1.1057)  =  0.5140  Btu/lbm  (15.48) 

e     p  s 

The  ideal  pump  work  is 
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*  s7    "   P6}  144       (1500.0    -    600.0) 

w 


p  P6  778.26  49.75 

=   3.3472  Btu/lbm  (15.49) 

The  net  useful  work  of  the  cycle  is  therefore 

*     * 
nxWe    Wp       (1) (0.5140)    3.3472 

wnet  '     x     n  0.4  72         0.9 

P 

=  -  2.6301   Btu/lbm  (15.50) 

Notice  that  the  net  work  of  the  cycle  happens  to  be  negative 
at  this  particular  operating  point.   In  other  words  the  electrical 
power  output  is  insufficient  to  drive  the  pump  under  these  condi- 
tions. 

The  net  heat  input  to  the  cycle  is 

* 

q..      =       (h   ,    -    hj E-  =    673.9    -    3'l4l2 

^■m  si  6'         n  0.9 

P 

=   670.2   Btu/lbm  (15.51) 

Finally,  the  overall  thermal  efficiency  of  the  system  at  this 
operating  point  is 

nc   =  -^^  =  "^n3.1  =  -  0.003924  (15.52) 

q.       670  . 2 

Although  the  performance  at  this  arbitrary  operating  conditions 
is  very  poor,  the  fact  remains  that  there  exists  some  (as  yet 
unknown)  value  of  the  parameter  vi/vo  =  Y   f°r  which  the  efficiency 
n    reaches  its  maximum  possible  value. 
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This  optimum  value  of   y  can  be  determined  by  further 
systematic  trial  calculations  similar  to  the  above  sample 
calculation. 

The  foregoing  calculations  have  established  the  values  of 
Ts  ,  Pg   and  M  at  all  stations.   The  corresponding  values  of 
T  ,  P  ,  p   and   s   at  these  locations  can  then  be  determined 
from  the  following  relations. 


T 

^  ■  1  +  (Y  ~  1)M2  =  1  +  0.1098M2  (15.53) 

T  =  Tg/(TS/T)  °R  (15.54) 


P_     T   y/(Y  -  1)     T   5.5537 

p^  =  (t1)  =  <^>  (15.55) 

P  =  P  /(P  /P)   psia  (15.56) 

P       1 44     p  "5 

P  =  ST  =  (F5TT5)  T     lbm/ft  (15.57) 

T  P 

s  =  Cln    (=2-)  -  RJln(=£)  (15.58) 

p     o         o 

T  P 

=  0.4649£n(5^Q  Q)  -  0 . 08371£n (^    696)  Btu/lbm  °R     (15.59) 

The  calculation  of  the  above  properties  has  been  completed 
in  this  way  and  the  results  are  summarized  in  Table  15.1.   These 
same  states  are  shown  on  a  temperature  entropy  diagram  in  Fig.  15.1 

Notice  from  the  tabulation  that  in  passing  through  the 

3 
ejector,  the  fluid  density   p   changes  from  1.402  lbm/ft   to 

3 
1.524  lbm/ft   ,  an  increase  of  about   8.5%  .   In  passing  through 

the  working  section,  the  density  changes  from   1.524  lbm/ft   to 


79 


w 

H 

Eh 
« 

cu 
o 

u 

H 

>h 
Q 
O 


W 

M 

En 
in 


CN 


MD 


CN 

• 

o 

in       m 
o 


ro 


m 
o 


ro 


m        ro 
o 


m 


>i      ro 

o 


ro 


ro 
o 


10 

O 

•<T 

O 

en 

\o 

« 


o 

rH 
CM 


ro 
oo 


in 
r- 


in 
r- 


m 

CTi 


en 


ro 

CO 


o 


o 


CN 


m 
o 

rH 

o 


en 


o 

rH 


CTi 


o 

rH 


Cn 

• 

CN 

H 

O 


m 

r* 

CN 

• 

(N 

en 

ro 

rH 

• 

o 

o 
o 


CO 

CM 
Cn 


en 

CTi 


Cn 

• 

cn 


in 
<-\ 


CO 

o 
en 


in 
e'- 


en 


in 

rH 

in 


CN 

m 


CM 

in 


m 


o 
r- 

CTi 


cu 


ro 


O 

I 


ro 


o 
O 


r- 
ro 


o 
I 


in 
m 


o 
o 


CN 
CN 

ro 
o 

o 

l 


in 

CN 


CN 


CTi 

H 

CN 


K 


g 
X! 


3 
Co- 
rn 

r- 
ro 

oo 
o 


PS 


g 

i2 


3 
CQ 
en 


I 


>i     X     >~    tf 


u 


a 


CO 

r> 

vo 

Cn 
II 
CO. 


w 

J 

r> 

CQ 

>£> 

r-» 

o 

CN 

O 

CN 

CN 

< 

• 

• 

o 

• 

• 

O 

CN 

Eh 

CN 

m 

ro 

ro 

<D 

o 

•*r 

ro 

m 

ro 

• 

■"^ 

o 

• 

O 

CTi 

U3 

o 

Cn 

V£> 

rH 

O 
1 

Cn 

O 

• 

• 

rH 

CN 

O 

CN 

in 

o 

00 

• 

• 

CN 

CN 

rH 

rH 

o 

r-~ 

^O 

o 

o 

ro 

rH 

in 

CN 

•^ 

o 

■^ 

o 

rH 

rH 

• 

rH 

Cn 

<o 

rH 

• 

o 
1 

a 


ro 

g 

4J 

Xi 

IH 

rH 

fd 

\ 

\ 

•H 

g 

3 

e< 

W 

JQ 

4-> 

o 

a 

rH 

CQ 

H 

CU 

Q. 

W 

80 


-i     o 


CO 

ro    co^j- 

co     ,    • 

^  ♦  ro 


>■ 

CO 


ro 

CO 


CO" 


M 

XX 

CO 


"A 

CVJ      CVJ 

co    ^-~ 


L 

O 
O 

C\J 


I 


I 


CO 

UJ 

CO 

o 


>- 

Q 

^^ 

O 

O 

5 

d 

Ul 

i 

X 

i- 

u. 
o 

— 

2 

cr 

o 

< 
cr 

< 

E 

-Q 

CVJ 

Q 

\   " 

o 
b 

>- 

a. 

3 

i 

o 

z 

a: 

z 

LU 

> 

a.  - 

UJ 

o 

o: 

or 

3 

I- 

}- 

z 

< 

LU 

01 

ro 

UJ 

O 

Q. 

d 

2 
UJ 

1 

1- 

m 
u. 


o 
o 


o 
o 
o 


o 
o 

CD 


o 
o 

00 


o 

o  d 

o    ■ 


Uo    3Un±VU3dW3± 


81 


3 

1.515  lbm/ft   ,  a  decrease  of  well  under   1%  .   These  data  support 

the  principle  that  it  is  permissible  to  neglect  the  density  changes 
across  these  components  for  purposes  of  simplifying  the  analysis. 
This  is  particularly  justified  because  at  present  the  uncertainity 
in  the  value  of  the  empirical  ejector  effectiveness  r\„      far 
exceeds  the  small  error  introduced  by  neglecting  these  minor  density 
variations.   Moreover,  the  present  sample  calculation  applies  to  a 
velocity  ratio  vi/vo   °^   4.25   which  is  an  extreme  case.   For 
efficient  ejector  operation  the  ratio  Vi/V9   should  be  only 
slightly  greater  than  unity  and  under  these  circumstances  the  den- 
sity change  across  the  ejector  is  expected  to  be  far  less  than 
8%  .   Thus  there  is  little  substance  in  the  argument  that  an 
analysis  which  neglects  these  minor  compressibility  effects  must 
be  seriously  in  error.   Such  an  argument  runs  counter  to  the 
fundamental  principles  of  fluid  mechanics.   A  more  correct  appraisal 
of  the  situation  is  that  the  inclusion  of  these  detailed  compressi- 
bility effects  can  be  expected  to  improve  the  accuracy  of  the 
analysis  to  some  degree,  but  that  they  should  not  be  expected  to 
change  the  overall  results  in  any  radical  way. 
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16.   Nomenclature 

A    -  area 

a    -  speed  of  sound 

C0    -  characteristic  breakdown  constant  (Eqn  9.3) 

a 

C    -  specific  heat  at  constant  pressure 

cf   -  friction  loss  coefficient  -  condenser 

E,    -  breakdown  electric  field  strength 

E(x)  -  Eqn.  10.13 

F(M  )-  Eqn.  4.16 

f(M)  -  Eqn.  4.3 

f     -  Eqn.  4.15 

max  ^ 

g(M)  -  Eqn.  6.5 

g     -  Eqn .  6.7 

q     -  Eqn.  6.9 

^max  ^ 


9oo 

- 

Eqn.  6.10 

h 

- 

enthalpy 

I 

- 

irreversibility 

Ko 

- 

Eqn.  10.11 

Kl 

- 

Eqn.  10.12 

M 

- 

Mach  number 

m 

- 

mass  flow  rate 

Pe 

- 

gross  electrical  power 

P 

- 

pressure 

qin 

- 

heat  input 

R 

- 

gas  constant 

R 

— 

universal  gas  constant 
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s  -  entropy 

s  -  entropy  of  mixing 

T  -  temperature 

V  -  velocity 

W  -  molecular  weight 

* 

w  -  gross  work  out  per  unit  mass 

* 

w  -  ideal  pump  work  per  unit  mass 

w  useful  electrical  output  per  unit  mass 
net  r   c 

x  -  itu/m- 

xnew  '  E<3"-    10"15 

y  -  v1/v3 

z  -  quality  of  vapor 

B  -  Eqn.  9.7,  dimensionless  breakdown  constant 

y  -  ratio  of  specific  heats 

A (  )  -  change  of 

e  -  permittivity  of  the  medium 

n  -  overall  cycle  efficiency 

t)  -  Eqn.  7.13,  ejector  effectiveness 
E 

nn  -  Eqn.  8.6,  diffuser  effectiveness 

D  -  pump  efficiency 

ri  -  excitation  efficiency 

p  -  gas  density 

iP  -  steady  flow  availability 

w  -  mass  ratio  of  condensible  primary  to  dry  secondary 
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Subscripts 

0  -  ambient  conditions 

1  -  primary,  entrance  to  ejector  (exit  from  boiler/superheater) 

2  -  secondary,  entrance  to  ejector  (exit  from  separator) 

3  -  mixture,  exit  from  ejector  entrance  to  diffuser 

4  -  mixture,  exit  from  diffuser  entrance  to  generator 

5  -  mixture,  exit  from  generator  entrance  to  condenser 

6  -  primary,  exit  from  condenser  entrance  to  pump 

7  -  primary,  exit  from  pump  entrance  to  boiler 
f  -  saturated  liquid 

G  -  gas  or  secondary 

g  -  saturated  vapor 

n  -  next  value 

s  -  stagnation 

v  -  vapor 

x  -  exit  station  of  hypothetical  minimum  loss  device 

y  -  exit  station  of  ideal  ejector 


85 


17.   References 

1.  T.  H.  Gawain  and  0.  Biblarz,  "Performance  Analysis  of  a  Type 
of  Electrohydrodynamic  Power  Generator,"  NPS67-79-006  (April 
1979) . 

2.  T.  H.  Gawain  and  0.  Biblarz,  "Performance  Analysis  of  a  Type 
of  Electrohydrodynamic  Generator,"  Proceedings  of  the  14th 
IECEC  (Boxton,  August  1979),  pp.  162-167. 

3.  T.  H.  Gawain  and  0.  Biblarz,  "Ejector-Augmented  Electro- 
hydrodynamic Power  Generator,"  submitted  to  J.  of  Chem.  Engr. 

4.  N.  M.  Huberman,  H.  Shelton,  W.  Krieve,  and  C.  L.  Dailey, 
"Study  of  Electrofluid  Dynamic  Power  Generation,"  AFAPL-TR- 
76-31  (July  1976) . 

5.  G.  B.  Wallis,  One-Dimensional  Two-Phase  Flow,  McGraw-Hill, 
N.  Y.  (1969) . 

6.  S.  H.  Hasinger,  "Comparison  of  Experiment  and  Analysis  for  a 
High  Primary  Mach  Number  Ejector,"  AFFDL-TR-38  (May  1977); 
"Ejector  Optimization,"  AFFDL-TR-23  (June  1978). 

7.  B.  H.  Anderson,  "Factors  which  Influence  the  Analysis  and 
Design  of  Ejector  Nozzles,"  AIAA  Paper  No.  72-46  (1972). 

8.  D.  E.  Groesbeck,  R.  G.  Huft,  and  U.  H.  von  Glahn,  "Comparison 
of  Jet  Mach  Number  Decay  Data  with  Correlation  and  Jet 
Spreading  Contours  for  a  Large  Variety  of  Nozzles",  NASA  TN 
D-8423  (1977) . 


86 


DISTRIBUTION  LIST 


No.  of  Copies 


1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  VA   22314 

2.  Library  2 
Code  0212 

Naval  Postgraduate  School 
Monterey,  CA   93940 

3.  Dean  of  Research  1 
Code  012 

Naval  Postgraduate  School 
Monterey,  CA   93940 

4 .  Chairman  1 
Department  of  Aeronautics 

Code  67 

Naval  Postgraduate  School 

Monterey,  CA   93940 

5.  Dr.  Ryszard  Gajewski,  Director  1 
Division  of  Advanced  Energy  Projects 

Mail  Stop  J-309 
Department  of  Energy 
Washington,  DC   20545 

6.  Mr.  Alvin  Marks,  Chief  Scientist  2 
Marks  Polarized  Corporation 

153-16  10th  Avenue 
Whitestone,  NY   11357 

7.  Professor  T.  H.  Gawain  2 
Department  of  Aeronautics 

Code  67Gn 

Naval  Postgraduate  School 

Monterey,  CA   93  940 

8.  Professor  0.  Biblarz  3 
Department  of  Aeronautics 

Code  67Bi 

Naval  Postgraduate  School 

Monterey,  CA   93940 


87 


No.  of  Copies 

9.   Dr.  H.  R.  Rosenwasser  1 

Naval  Air  Systems  Command 
Code  AIR  310C 
Washington,  D.  C.   203  60 

10.  Mr.  S.  A.  Satkowski  1 
Office  of  Naval  Research 

Power  Program,  Code  473 
Washington,  D.  C.   20360 

11.  Dr.  S.  Hasinger  1 
Thermomechanics  Branch 

AFFDL 

Wright-Patterson  AFB,  Ohio   45433 

12.  Dr.  C.  D.  Hendricks  1 
Lawrence  Livermore  Laboratory 

L-482 

P.  0.  Box  808 

Livermore,  CA   94550 

13.  Dr.  H.  Velkoff  1 
Department  of  Mechanical  Engineering 

Ohio  State  University 
Columbus,  Ohio   43210 

15.  Dr.  Hans  von  Ohain  1 
c/o 

University  of  Dayton  Research  Institute 
Dayton,  Ohio   45469 

16.  Dr.  J.  R.  Melcher  , 
Department  of  Electrical  Engineering 

Room  36-319 

Massachusetts  Institute  of  Technology 

Cambridge,  MA   02139 

17.  Dr.  J.  E.  Minardi  1 
University  of  Dayton  Research  Institute 

Dayton,  Ohio   45409 


88 


DUDLEY  KNOX  LIBRARY  -  RESEARCH  REPORTS 

5  6853  01068116  6 


